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Abstract. We develop a renormalisation scheme for time-ordered products in interacting field 
theories on curved spacetimes which consists of an analytic regularisation of Feynman amplitudes and 
a minimal subtraction of the resulting pole parts. This scheme is directly applicable to spacetimes 
with Lorentzian signature, manifestly generally covariant, invariant under any spacetime isometries 
present and constructed to all orders in perturbation theory. Moreover, the scheme captures correctly 
the non-geometric state-dependent contribution of Feynman amplitudes and it is well-suited for 
practical computations. To illustrate this last point, we compute explicit examples on a generic curved 
spacetime, and demonstrate how momentum space computations in cosmological spacetimes can be 
performed in our scheme. In this work, we discuss only scalar fields in four spacetime dimensions, but 
we argue that the renormalisation scheme can be directly generalised to other spacetime dimensions 
and field theories with higher spin, as well as to theories with local gauge invariance. 
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1. Introduction 

In the perturbative construction of models in quantum field theory on curved spacetimes one encounters 
time-ordered products of field polynomials which are a priori ill-defined due to the appearance of UV 
divergences. Several renormalisation schemes which deal with these divergences in the presence of non¬ 
trivial spacetime curvature have been discussed in the literature, such as for example local momentum 
space methods [Bu81], dimensional regularisation in combination with heat kernel techniques [Lii82, 
To82], differential renormalisation [CHL95, Pr97], zeta-function renormalisation [BF13], generic 
Epstein-Glaser renormalisation [BFOO, HW01, HW04], and, on cosmological spacetimes, Mellin-Barnes 
techniques [HolO] and dimensional regularisation with respect to the comoving spatial coordinates 
[BCK10], 

Some of these schemes, such as heat kernel approaches, zeta-function techniques and local momen¬ 
tum space methods are based on constructions which are initially only well-defined for spacetimes with 
Euclidean signature. These constructions can be partly transported to general Lorentzian spacetimes 
by local Wick-rotation techniques developed in [Mo99]. However, whereas the Feynman propagator 
is essentially unique on Euclidean spacetimes, this is not the case on Lorentzian spacetimes where 
this propagator has a non-unique contribution depending on the quantum state of the field model. 
Consequently, the Euclidean renormalisation techniques, and the numerous practical computations 
already performed by means of these methods - see for example the monographs [BD82, BOS92, PT09] 
- are able to capture the correct divergent and geometric parts of Feynman amplitudes, but a priori 
not their non-geometric and state-dependent contributions. 

A renormalisation scheme which is directly applicable to curved spacetimes with Lorentzian 
signature has been developed in [BFOO, HW01, HW04] in the framework of algebraic quantum field 
theory. This scheme implements ideas of [EG73] and [St71] and is based on microlocal techniques which 
replace the momentum space methods available in Minkowski spacetime and have been introduced 
to quantum field theory in curved spacetime by the seminal work [Ra96]. However, although the 
generalised Epstein-Glaser scheme developed in [BFOO, HW01, HW04] is conceptually clear and 
mathematically rigorous, it is not easily applicable in practical computations. On the other hand, 
Lorentzian schemes which are better suited for this purpose have not been developed to all orders in 
perturbation theory [CHL95, Pr97], are tailored to specific spacetimes [HolO] or are not manifestly 
covariant [BCK10]. 

Motivated by this, we develop a renormalisation scheme for time-ordered products in interacting 
field theories on curved spacetimes which is directly applicable to spacetimes with Lorentzian signature, 
manifestly generally covariant, invariant under any spacetime isometries present and constructed 
to all orders in perturbation theory. Moreover, the scheme captures correctly the non-geometric 
state-dependent contribution of Feynman amplitudes and it is well-suited for practical computations. 
In this work, we discuss only scalar fields in four spacetime dimensions, but we shall argue that the 
renormalisation scheme can be directly generalised to other spacetime dimensions and field theories 
with higher spin, as well as to theories with local gauge invariance. Our analysis will take place in the 
framework of perturbative algebraic quantum field theory (pAQFT) [BFOO, HW01, HW04, BDF09, 
FR12, FR14] which is a conceptually clear framework in which fundamental physical properties of 
perturbative interacting models on curved spacetimes can be discussed. However, we will make an effort 
to review how the formulation of pAQFT is related to the more standard formulation of perturbative 
QFT. 

The renormalisation scheme we propose is inspired by the works [KelO, DFKR14] which deal 
with perturbative QFT in Minkowski spacetime. In these works, the authors introduce an analytic 
regularisation of the position-space Feynman propagator in Minkowski spacetime which is similar to 
the one discussed in [BG72]. Based on this, time-ordered products are constructed recursively by an 
Epstein-Glaser type procedure and it is shown that this recursion can be resolved by a position-space 
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forest formula similar to the one of Zimmermann used in BPHZ renormalisation in momentum space. 

In order to extend the scheme proposed in [DFKR14] to curved spacetimes, and motivated by 
[BG72] and by the form of Feynman propagators on curved spacetimes, we introduce an analytic 
regularisation of a Feynman propagator Ap by 

:= lim (j -- + - (l- 7 - Ky-Y) + u>, 

e->o+87 t 2 \(cr + *e) 1+ “ a \ (a + ie) a J J 

where u, v and w are the so-called Hadamard coefficients and a is 1/2 times the squared geodesic 
distance. This analytic regularisation, namely the construction of certain distributions by means of 
powers of the geodesic distance, is reminiscent of the use of Riesz distributions to define advanced and 
retarded Greens functions on Minkowski spacetime. A careful discussion of Riesz distributions and 
their extension to the curved case is presented in [BGP07]. The regularisation we use is loosely related 
to dimensional regularisation because the leading singularity of a Feynman propagator in N spacetime 
dimensions is proportional to (a + ie) 1 ^ N ^ 2 1 see e.g. [Mo03, Appendix A]. A regularisation of the 
Feynman propagator similar to the one above has recently been discussed in [Dal5]. In this work, we 
shall combine the analytic regularisation of the Feynman propagator with the minimal subtraction 
scheme encoded in a forest formula of the kind discussed in [HolO, KelO, DFKR14] in order to obtain 
a time-ordered product which satisfies the causal factorisation property, i.e. a product which is indeed 
“time-ordered”. In order to prove that the analytically regularised amplitudes constructed out of A ^ 
have the meromorphic structure necessary for the application of the forest formula and in order to show 
how the corresponding Laurent series can be computed explicitly, we shall make use of generalised 
Euler operators. The practical feasibility of the renormalisation scheme shall be demonstrated by 
computing a few examples. 

A large part of our analysis is devoted to demonstrating that the scheme we propose is consistent 
and well-defined on general curved spacetimes. Readers interested in directly applying our scheme 
may use Proposition 3.9 with (3.22) and (3.23) in order to compute the Laurent-series of the Feynman 
amplitudes (3.11) regularised by means of (3.12). The correct order of pole subtractions is encoded in 
the forest formula (3.9) which is explained prior to its display. A number of examples is discussed in 
Section 3.5. 

In quantum held theory on cosmological spacetimes, i.e. Friedmann-Lemaitre-Robertson-Walker 
(FLRW) spacetimes, one usually exploits the high symmetry of these spacetimes in order to evaluate 
analytical expressions in spatial Fourier space. However, the renormalisation scheme discussed in 
this work operates on quantities such as the geodesic distance and the Hadamard coefficients, whose 
explicit position space and momentum space forms are not even explicitly known in FLRW spacetimes. 
Notwithstanding, we shall devote a large part of this work in order to develop simple methods to 
evaluate quantities renormalised in our scheme on FLRW spacetimes in momentum space, and we 
shall illustrate these methods by explicit examples. 

The paper is organised as follows. In the next section we present a brief introduction to pAQFT 
and its connection with the more standard formulation of perturbative QFT. Afterwards we introduce 
the renormalisation scheme, demonstrate that it is well-defined and analyse its properties in Section 3, 
where we also illustrate the scheme by computing examples. In the fourth section we demonstrate 
the applicability of the renormalisation scheme to momentum space computations on cosmological 
spacetimes. Finally, a few conclusions are drawn in the last section of this paper. Conventions regarding 
the various propagators of a scalar field theory and a few technical computations are collected in the 
appendix. 


2. Introduction to pAQFT 

2.1. Basic definitions 

Throughout this work, we shall consider four-dimensional globally hyperbolic spactimes (A4,g), where 
g is a Lorentzian metric whose signature is (—,+,+,+) and we use the sign conventions of [Wa84] 
regarding the definitions of curvature tensors. 


3 




We recall the perturbative construction of an interacting quantum field theory on a generic curved 
spacetime in the framework of perturbative algebraic quantum field theory (pAQFT) recently 
developed in [BDF09, FR12, FR14] based on earlier work. In this construction, the basic object 
of the theory is an algebra of observables which is realised as a suitable set of functionals on field 
configurations equipped with a suitable product. In order to implement the perturbative constructions 
following the ideas of Bogoliubov and others, the field configurations 4> are assumed to be off-shell. 
Namely, 4> £ £{M.) = G°°( M) is a smooth function on the globally hyperbolic spacetime (A i,g) and 
observables are modelled by functionals F : £(Ad) C satisfying further properties. In particular all 
the functional derivatives exist as distributions of compact support, where we recall that the functional 
derivative of a functional F is defined for all Vd, ■ • •, V’n £ V(M.) = as 

fjn 

F^(mi ®---®ipn):= 77 w F{<j> + ArV-i + • • • A„V>n) £ £'(M n ). 

aAi ... a\ n Ai=—=A„=o 

The set of these functionals is indicated by F . Further regularity properties are assumed for the 
construction of an algebraic product. In particular, the set of local functionals F\ oc C T is formed 
by the functionals whose n—th order functional derivatives are supported on the total diagonal 
d n = {(x,..., x), x £ A4} C Ad”. Furthermore, their singular directions are required to be orthogonal 
to d n , namely WF(F(”)) C {(a:, k) £ T*M n ,x € d n ,k _L Td n } where WF denotes the wave front 
set. A generic local functional is a polynomial P{4>){x) in <f> and its derivatives integrated against a 
smooth and compactly supported tensor. The functionals whose functional derivatives are compactly 
supported smooth functions are instead called regular functionals and indicated by F reg . 

The quantum theory is specified once a product among elements of J-i oc and a *—operation (an 
involution on F) are given. For the case of free (linear) theories the product can be explicitly given by 
a ^product 

Un . . 

F* H G=Y J — (F {n \Hf n G( n) \ , ( 2 . 1 ) 

n 

where H + is a Hadamard distribution of the linear theory we are going to quantize, namely a 
distribution whose antisymmetric part is proportional to the commutator function A = Ar — A a and 
whose wave front set satisfies the Hadamard condition, see e.g. [Ra96, BFK96] for further details and 
Section A.l for our propagator conventions. Owing to the properties of H +1 iterated *y-products of 
local functionals F\ *h • • • *h F n are well defined and is associative. 

In a normal neighbourhood of (A4 , 3 ), a Hadamard distribution H + is of the form 

H+(x,y) = ( U ^ ,yS) + v{x,y)log{M 2 a + (x,y))] +w{x,y), (2.2) 

8 tH \a+{x,y) J 

where cr+(x,y) = cr(x,y) + ie(t(x) — t{y)) + e 2 /2 with t a time-function, i.e. a global time-coordinate, 
2a(x,y) is the squared geodesic distance between x and y and M is an arbitrary mass scale. The 
Hadamard coefficients u and v are purely geometric and thus state-independent, whereas w is smooth 
and state-dependent if H + (x,y) is the two-point function of a quantum state. 

For the perturbative construction of interacting models we further need a time—ordered product 
■t h on local functionals. This product is characterised by symmetry and the causal factorisation 
property, which requires that 


F - Th G = F -k H G if F > G, (2.3) 

where F > G indicates that F is later than G, i.e. there exists a Cauchy surface E of (A 4,g) such that 
supp(F) C T + (E) and supp(G) C J“(E). However, the causal factorisation fixes uniquely only the 
time-ordered products among regular functionals, in which case 

F - Th G — Y^^{ FM ,Hf n G^ , (2.4) 

n 

where Hr is the time-ordered (Feynman) version of H + , i.e. Hr = H + + iA a with A^ the advanced 
propagator of the free theory, cf. Section A.l. For local functionals, (2.4) is only correct up to the 
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need to employ a non-unique renormalisation procedure, cf. Section 3.1. This renormalisation can 
be performed in such a way that iterated -t h -products of local functionals F\ - Th ■ ■ ■ - Th F n are well 
defined with -t h being associative. Moreover, -products of such time-ordered products of local 
functionals are well-defined as well, cf. [HW02, BDF09, FR12, FR14]. Consequently, we may consider 
the algebra Ao *//-generated by iterated -^-products of local functionals. This algebra contains all 
observables of the free theory which are relevant for perturbation theory. 

In the perturbative construction of interacting models, namely when the free action is perturbed by 
a non-linear local functional V, the observables associated with the interacting theory are represented 
on the free algebra Ao by means of the Bogoliubov formula. This is given in terms of the local 
S'-matrix, i.e., the time-ordered exponential 

oo . n 

= P-5) 

n =0 

n times 

where V is the interacting Lagrangean. In particular, for every interacting observable F the corre¬ 
sponding representation on the free algebra Aq is given by 

7 lv(F) = S-\V ) ( S(V ) -Th F) , (2.6) 

where S ,_1 (V r ) is the inverse of S(V) with respect to the *ff-product. The problem in using TZy(F) 
as generators of the algebra of interacting observables lies in the construction of the time-ordered 
product which a priori is an ill-defined operation. 

This problem can be solved using ideas which go back to Epstein and Glaser, see e.g. [BFOO], by 
means of which the time-ordered product among local functionals is constructed recursively. The 
time-ordered products can be expanded in terms of distributions smeared with compactly supported 
smooth functions which play the role of coupling constants (multiplied by a spacetime-cutoff). At 
each recursion step the causal factorisation property (2.3) permits to construct the distributions 
defining the time-ordered product up to the total diagonal. The extension to the total diagonal can be 
performed extending the distributions previously obtained without altering the scaling degree towards 
the diagonal. In this procedure there is the freedom of adding finite local contributions supported 
on the total diagonal. This freedom is the well known renormalisation freedom. In addition to the 
properties already discussed, the renormalised time-ordered product is required to satisfy further 
physically reasonable conditions. We refer to [HW02, HW04] for details on these properties and the 
proof that they can be implemented in the recursive Epstein-Glaser construction. 

In spite of the theoretical clarity of this construction, the Epstein-Glaser renormalisation is quite 
difficult to implement in practise. The aim of this paper is to discuss a renormalisation scheme which 
is suitable for practical computations. 

2.2. Relation to the standard formulation of perturbative QFT 

In this subsection we outline the relation of the pAQFT framework to the standard formulation of 
perturbative QFT. As an example, we demonstrate how the two-point (Wightman) function of the 
interacting field in cjA theory on a four-dimensional curved spacetime is computed, where we assume 
that the quantum state of the interacting field is just the state of free field modified by the interacting 
dynamics. We further assume that the free field is in a pure and Gaussian Hadamard state. 

Let us recall the relevant formulae in perturbative algebraic quantum field theory where we shall 
always try to write expressions both in the pAQFT and in the more standard notation, indicating the 
latter by a =. Given a local action V, such as V = f M d 4 x-y^g j4>{x) 4 in </> 4 -theory, the corresponding 
S-matrix, which is loosely speaking the “S'-matrix in the interaction picture”, is defined by (2.5) and 
corresponds to S(V) =Te^ v . 

The interacting field, i.e. “the field in the interaction picture” (f>i{x), is defined by the Bogoliubov 
formula 

M*) = Rv^(i)) = SiV)- 1 * H ( S(V ) • Th <Kx)) = T{e% v )~ 1 T{e^ v (2.7) 

similarly to (2.6), where by unitarity S'(E) -1 = S(V)*. Interacting versions of more complicated 
expressions in the field, e.g. polynomials at different and coinciding points, are defined analogously. A 
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thorough discussion of the relation between the Bogoliubov formula and the more common formulation 
of observables in the interaction picture may be found e.g. in [Lil3, Section 3.1]. We only remark that, 
in the Minkowski vacuum state the expectation value of the Bogoliubov formula can be shown to 
read (also for more general expressions in the held) 


(M*))n 0 = (ne^r^ei^Or)))^ 


\ / S 21 


which is the theorem of Gell-Mann and Low, see [Du96, DFOO] for details. 

In the algebraic formulation one usually cuts off the interaction in order to avoid infrared problems 
by replacing A — > Xf(x) with a compactly supported smooth function / and considers the adiabatic 
limit / —> 1 in the end when computing expectation values. As our aim is to compute expectation 
values in this section, we shall write the results in the adiabatic limit keeping in mind that proving the 
absence of infrared problems, i.e. the convergence of the spacetime integrals, is non-trivial and may 
depend on the state of the free held chosen. Note that the so-called “in-in-formalism” often used in 
perturbative QFT on cosmological spacetimes corresponds to considering a cutoff function / of the 
form f(t, x) = @(f — to), i.e. / is a step function in time and the parameter to corresponds to the time 
where the interaction is switched on. 

Our choice for the quantum state O of the interacting held implies that e.g. the interacting 
two-point function 

(&i(x) * H <t>i{y))n = {(/) I (x)(j) I (y))n 


is computed by writing </>/ in terms of the free held <fi and computing the expectation value of the 
resulting observable of the free held in the chosen pure, Gaussian, Hadamard state of the free held which 
we may thus denote by the same symbol fi. The interacting vacuum state in Minkowski spacetime is 
of this form, whereas interacting thermal states in hat spacetime do not belong to this class, as they 
roughly speaking require to take into account both the change of dynamics and the change of spectral 
properties induced by V [FL13]. 

The functionals in the functional picture of pAQFT correspond to Wick-ordered quantities of 
the free held in the sense we shall explain now. To this avail we recall the form of the (quantum) 
*h product and (time-ordered) -^-product in (2.1) and (2.4) which are dehned by means of a 
Hadamard distribution H + and its Feynman-version H F = H + + iA a- Up to renormalisation of the 
time-ordered product, these products computed for the special case of the functional <fi 2 (x) give 


</>(z) 2 * H </>(y) 2 = fiix) 2 4>{y) 2 + 4 h(j){x)(j){y)H + {x, y) + 2 H 2 Hl(x, y ), 

<t>(x) 2 - Th </>(y) 2 = <Ka;) 2 </>(y) 2 + 4. hcj)(x)<t>{y)H F {x , y ) + 2h 2 H 2 F (x, y). 

This example shows that the product ( -t h product) implements the Wick theorem for normal- 
ordered (time-ordered) helds, and thus the previous formulae can be interpreted in more standard 
notation as 


: (j){x) 2 : H : (j){y) 2 : H =: 4>(x) 2 (j){y) 2 : H +4 h : <^( x)4>{y) : H H + (x, y) + 2Ti 2 # 2 (z, y) , 

T (: 4>(x) 2 : H ■ <t>{y) 2 -h) =■ c/)(x) 2 cj)(y) 2 : H +4fi : <!>(x)(j>(y) : H H F (x, y) + 2h 2 H F {x, y ), 

where 


:A:h := a-H + (A) := e n ( Hs ( x ’ v ' > ’/)) a, 
H s (x, y) := i( H+{x , y) + H + {y, x)), 

e.g. 


( 2 . 8 ) 


:<p(x) 2 : H = lim (<j>{x)^{y) - H+(x,y)) . 

x^y 

(2.8) is a convenient way to encode the combinatorics of normal ordering whereby the exponential 
series terminates for polynomial functionals such as A = (f>(x) 2 . 
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The Wick theorem relates (time-ordered) products of Wick-ordered quantities to sums of Wick- 
ordered versions of contracted products, where the definition of “Wick-ordering” and “contraction” 
are directly related, they both depend on the Hadamard distribution H + chosen. Thus, if we choose a 
particular H + to define *h and -t h in pAQFT, we immediately fix the interpretation of all functionals 
in terms of expressions Wick-ordered with respect to H + . 

For the algebraic formulation the choice of H + is not important, indeed choosing a different H' + 
with the same properties, one has that w := H' + — H + = H' F — Hp , because the advanced propagator 
A 4 is unique and thus universal. Moreover, w is real, smooth and symmetric and 

A'kH' dB — cx w (o_u,(A) -kp , A 't h i dB — a w (^—^(A) ‘t h w(dB ')'), 

with a defined as in ( 2 . 8 ) and thus the algebras associated to *#, -t h and -t h , are isomorphic via 

a w '■ Aq —> Aq , 

where we recall that Aq is algebra *#-generated by -^-products of local functionals. 

Flence, one may choose a suitable H + according to ones needs. Flowever, since ad(A) ^ A for 
functionals containing multiple field powers, statements like “the potential is </> 4 ” are ambiguous in 
pAQFT, and in fact also in the standard treatment of QFT. They become non ambiguous only if one 
says “the potential is :</> 4 :#, i.e. (f> 4 Wick-ordered with respect to H + ”. In pAQFT the corresponding 
non-ambiguous statement would be “the potential is the functional (jA in the algebra Ao constructed 
by means of dd+”. If one then passes to the algebra Aq constructed by means of H' + , the potential picks 
up quadratic and c-number terms as we shall compute explicitly below. Alternatively, this ambiguity 
may be seen to correspond to the renormalisation ambiguity of tadpoles in Feynman diagrams. 

Given a Gaussian and Hadamard free field state $7, a convenient choice or representation of the 
algebra is to take H + = A + , where A + (x, y) = (4>(x)*a </>(j/))o = (<Kx)<?i(?/))n is the two-point function 
of the free field in the state Cl. This corresponds to standard normal-ordering and consequently in this 
representation the expectation values of all expressions which contain non-trivial powers of the field 
vanish, i.e. 

(A)n = A\ ( f >=0 = (: A:a)o • (2.9) 

Keeping the state 12 fixed, but passing on to a representation of the algebra with arbitrary H + , the 
expectation value is computed as 

(A)n = aw(A) 1^=0 = (: A: H )n , w = A + — H + , 


for instance 


(^ 2 (aO)n = ct w (4> 2 (x ))|^=o = ( 0 2 (x) + w(x,x)) |^, =0 = w(x,x) = (:^ 2 (x) : H )n , 
which in more standard terms would be computed as 

(:cp 2 {x): H )n = lim {<t>{x)<t>{y) - H+(x, y)) n = lim (A +(x,y) - H+(x,y)) =w(x,x). 

x—>y x—>y 

In QFT in curved spacetimes normal-ordering is in principle problematic, because (pointlike) 
observables should be defined in a local and generally covariant way, i.e. they should only depend on 
the spacetime in an arbitrarily small neighbourhood of the observable localisation [BFV01, HW01]. 
This is not satisfied for e.g. field polynomials Wick-ordered with A + (x,y), because this distribution 
satisfies the Klein-Gordon equation and thus it encodes non-local information on the curved spacetime 
[HW01]. It is still possible to compute in the convenient normal-ordered representation in the following 
way. In the example of 0 4 -theory, one defines the potential j<f>(x) 4 as a local and covariant observable 
by identifying it with the corresponding monomial in a representation of the algebra furnished by a 
purely geometric H + , i.e. a H + of the form (2.2) with w = 0. 

In other words, we set once and for all in the .£/+-representation 

V H = [ d 4 x^g j^QO 4 == [ d 4 xy/^g j :c/)(x) 4 :h ■ 

Jm 4 Jm 4 
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This does not fix V uniquely, because H depends on the scale M inside of the logarithm, but the 
freedom in defining Vp, and analogously the free/quadratic part of Klein-Gordon action, as above 
corresponds to the usual freedom in choosing the “bare mass” to, “bare coupling to the scalar curvature” 
£, “bare cosmological constant” A, “bare Newton constant” G, as well as the “bare coefficients” pi, P 2 
of higher-derivative gravitational terms in the extended Einstein-Hilbert-Klein-Gordon action 


${<!>, 9ab)= [ d 

J M 




(R- 2A 
V 16t rG 


+ f3\ R 2 + /?2 RabR a — 


(v</> 2 ) 


(to 2 + £,R) 4> 2 

2 



In order to switch to the normal-ordered representation, we use the map a w defined in (2.8) where 
w = A + — H + is the state-dependent part of the Hadamard distribution A + whose dependence on the 
choice of M in H + corresponds to the above-mentioned freedom in the definition of the Wick-ordered 
Klein-Gordon action. That is, we have in the normal-ordered representation in the state f l 


r _ \ o\ o\ 

V :=V A = a w (V H ) = / d^Xy/^g 7 </>(x ) 4 + —w(x, x)(/)(x) 2 + —w(x, x) 2 (2.10) 

Jm 4 z 4 

= [ y :^(i) 4 : A + 7 ui(i,i) :(i(i) 2 : a + j«)(i,i ) 2 

Jm 4 z 4 

We observe that the combination of the requirements that the interaction potential is a local and 
covariant observable and that, in order to compute expectation values in the state 12 , one would like to 
compute in the convenient normal-ordered representation with respect to 12, leads to the introduction 
of an effective spacetime-dependent and state-dependent (squared) mass term g{x) = 3A w(x,x) in the 
interaction potential which of course leads to additional Feynman graphs in perturbation theory, cf. 
Figures 1 and 2. The field-independent term ^w(x, x) 2 plays no role for computations of quantities 
which do not involve functional derivatives of the extended Einstein-Hilbert-Klein-Gordon action 
with respect to the metric (an example where it does play a role is the stress-energy tensor), just as 
the modification of the free action by the change of representation plays no role for the computation 
of such quantities. A similar phenomenon as in (2.10) occurs in thermal quantum field theory on 
Minkowski spacetime, where the effective mass generated by changing from the normal-ordered picture 
with respect to the free vacuum state to the normal-ordered picture with respect to the free thermal 
state is termed “thermal mass”, cf. [Lil3, Section 2.3.2.] for details. 

After these general considerations, we can proceed to compute as an example the two-point function 
of the interacting field </>/ in </> 4 up to second order in A, whereby </>/ is assumed to be in a state induced 
by a Gaussian Hadamard state of the free field. To this avail, we shall exclusively compute in the 
associated normal-ordered representation and thus omit the subscripts on the star product, and the 
time-ordered product, * := *a, -t ■= t a - 

We start from the Bogoliubov formula (2.7) and compute (from now on h = 1) 

S(V) = 1 + iV - - T V + 0( A 3 ) 


S-(H )*” 1 = 1 -iV+^V- T V-V*V + 0( A 3 ) 

c/>I = (f> — iV -k (j) + iV -T 4> + - (V -T V) * <j) — V * V * (j) — -V -T V -T <t> + V * (V -t <j>) + 0( A 3 ). 

It remains to compute the *-product of (j>i{x) and 4>i{y) and to set (j) = 0 in the remaining expression 
in order to obtain the expectation value in the state f2. The result can as always be conveniently 
expressed in terms of Feynman diagrams, where we use the Feynman rules depicted in Figure 1. 

In the computation of (cf>i(x)(t>j(y)) o, many expressions can be shortened considerably by using 
the relation Ap ~ A + = iAa, in particular this holds for the external legs of the appearing Feynman 
diagrams. The resulting Feynman diagrams are depicted in Figure 2. 
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Figure 1. The various propagators and vertices in </> 4 -theory, where y(x) = 3A w(x,x). 
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Figure 2. The up-to-second-order contributions to the two-point (Wightman) function {<j>i(x)<j)i(y))n 
of the interacting field with potential j(f>(x) 4 + ^ —(f>(x ) 2 . We omit the labels of the external vertices 
after the first line using the convention that the left external vertex is always the x-vertex. 


3. Analytic regularisation and minimal subtraction on curved spacetimes 


As discussed above, the main problem in using the Bogoliubov formula (2.6) 


K V (F) 


-^S(V)~ X * H S(V + XF) 
i aX 


for constructing interacting fields perturbatively is that it is given in terms of the S'-matrix, which is 
the time-ordered exponential (2.5). Unfortunately, the time-ordered product defined in terms of a 
“deformation” (2.4) written by means of a Feynman propagator Hp is well defined only on regular 
functionals because the singularities present in Hp forbid their application to more general functionals. 

In order to proceed there is the need of employing a renormalisation procedure to construct the 
time-ordered products. In this work we discuss the use of certain analytic methods to solve this 
problem. The procedure we shall pursue is the following. We deform the Feynman propagator by 
means of complex parameter a with values in the neighbourhood of the origin obtaining a function 
with distributional values a i —> H^' 1 . The deformation we are looking for needs to be such that in the 
limit a —> 0 we recover the ordinary Feynman propagator. Furthermore, when a is non-vanishing, but 
sufficiently small, pointwise powers of Hp^ and integral kernels of more complicated loop diagrams 
should be well-defined. If this is the case, since the corresponding distributions obtained in the limit 
a —> 0 are well defined outside of the total diagonal, the poles of a <—> and more complicated loop 

expressions are supported on the total diagonal. The idea, similar to what happens in dimensional 
regularisation, is that it is possible to renormalise these distributions by simply removing the poles. 
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3.1. Analytic regularistion of time—ordered products and the minimal subtraction scheme 

In order to discuss the analytic regularisation of time-ordered products, we employ the notation used 
e.g. in [DFKR14] which efficiently encodes the full combinatorics of Feynman diagrams in a compact 
form. Namely, the time-ordered product of n local functionals V ±,..., V n can be formally defined in 
the following way 1 * 


V\ -t h ■ ■ ■ -t h Vn '■= T n (Vi ® ■ ■ ■ <8> V n ) := to o T n (Vi <8> ■ ■ ■ <8> V n ), (3.1) 

where m denotes the pointwise product m(F\ (g) • • • ® F n ){<j)) = F\(cj))... F n (<j>) and the operator T n is 
written in terms of an exponential 


( \ °° A hj 

T n = exp I £ Ay I = I] E]T ( 3 - 2 ) 

with 

A « := ( Hr ’ sm )' (33) 

Here the functional derivative acts on the i —th element of the tensor product V\ (g> • • • <g> V n 
and H F = H + + iAa is the time-ordered version of the Hadamard distribution Fl + entering the 
construction of the free algebra Ao via *#. The exponential (3.2) admits the usual representation in 
terms of Feynman graphs. More precisely, it can be written as a sum over all graphs T in Q n , the set of 
all graphs with vertices V(r) = {1,... ,n} and kj edges e € -E’(r) joining the vertices i,j. Furthermore, 
in this construction, there are no tadpoles h n = 0 (cf. Section 2.2 for details on why these are absent) 
and the edges are not oriented Uj = Iji. With this in mind 

rp ST' _j_ / 

r^g n ^ riiev(r) riB(r)9eDi 

where iV(r) = a numerical factor counting the possible permutations among the lines 

joining the same two vertices, the second product Y\ e0i is over the edges having * as a vertex and Xi 
is a point in M. corresponding to the vertex i. Moreover, ip is a distribution which is well-defined 
outside of all partial diagonals, namely on M n \ D n , where 

D n := {a;i,..., x n | Xi = Xj for at least one pair i ^ j} (3-5) 

and rr has the form 

Tr = H F (xi,Xj)= H F (xi,Xj) lij . (3.6) 

e=(i,j)EE( r) l<i<j<n 

The a priori restricted domain of rp is the reason why T n defined as above is not a well-defined 
operation on J 7 ®". In this context we recall that the total diagonal d n C D n is defined as 

d n := {(x,..., x), x € M} C M n . (3.7) 

In order to complete the construction we need to extend the obtained distributions to the diagonals 
D n . This is not a straightforward limit because the singular structure of the Feynman propagator Flp 
contains the one of the ^-distribution and because pointwise products of the latter distribution are 
ill-defined. Consequently, a renormalisation procedure needs to be implemented in order to extend rr 
to the full M. n . This extension is in general not unique, but subject to renormalisation freedom. 

1 In fact, in view of locality and covariance a better definition of the time-ordered product is 7i(Vi )-t h ■ ■ "T H Fi(V n ) '■= 

Tn(ViiS>- ■ -®V„) where 7i : T] oc —> T\ oc C Ao plays the role of identifying local and covariant (smeared) Wick polynomials 

as particular elements of the free algebra Ao, cf. [HW04]. As we shall not touch upon this point in our renormalisation 

scheme, we choose to omit 7i in our formulas for simplicity. 
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Here we shall discuss a procedure to extend the distributions rr to D n called minimal subtraction 
(MS), which makes use of an analytic regularisation A“ ,J of given in terms of a family of 
deformations Hp' J of the Feynman propagator Hp parametrised by complex parameters ctij contained 
in some neighbourhood of 0 £ C. To this end, we follow [DFKR14] and call tT) an analytic 
regularisation of a distribution t defined outside of a point x o £ M. if for all / £ V(A4) (t^ a \ f) is a 
meromorphic function in a for a in some neighbourhood of 0 which is analytic for a^O. Moreover 
t may be extended to Xq for a 0 whereas lim^^o t= t on Ad \ {cco}. 

We shall introduce an analytic regularisation of the Feynman propagator Hp in the following 
section, but the basic idea of the MS-scheme is independent of the details of the analytic regularisation. 

(ck) 

Namely, given any analytic regularisation Hp ' of Hp , we repeat the formal construction of T n 
presented above by replacing Hp by H^ in (3.3) and A^ by the induced A“ 1J in (3.2). Proceeding 
in this way we define 

:= e^i<j Ai i with a := {atijj^j , 

and the corresponding integral kernels r jT' 1 of Feynman graphs T in analogy to (3.4). We expect that 
the distributions are multivariate meromorphic functions which have poles at the origin for some 
of the ctij. Hence, in order to obtain well-defined distributions in the limit ctij to 0 and consequently 
a renormalised time-ordered product -t h , all these poles need to be subtracted. 

The properties of the analytically regularised Feynman propagator imply that T^ a ' > is well-defined 
on \ D n (3.5) even if all a l3 are vanishing. Since ris a multivariate meromorphic function 
in a which is analytic if restricted to M n \ D n , we may deduce that the principal part of for 
some aij must be supported on a partial diagonal of A4". In fact, in order for the time-ordered 
products to fulfil the factorisation property (2.3), the subtraction of the principal parts of needs 
to be done in such a way that at each step only local terms are subtracted. However, the previous 
discussion only implies that the support of the principal parts is contained in D n , i.e. the union of 
all the partial diagonals in AT 1 . In order to satisfy the causal factorisation property, the principal 
parts need to be removed in a recursive way starting from the partial diagonals corresponding to two 
vertices and proceeding with the partial diagonals corresponding to an increasing number m < n of 
vertices 0/ := {(aq,..., x n ) £ M n , aq = Xj,i,j £ I C {1,..., n}, |/| = m}. 

The correct recursion procedure is implemented by the so called Epstein-Glaser forest formula, 
which is a position-space analogue of the Zimmermann forest formula, see [HolO, KelO, DFKR14] for 
a careful analysis of the subject. We shall here follow the treatment discussed in [DFKR14], To this 
end, we consider the set of indices n := {1, ..., n} and define a forest F as 

F = {/i,...,4}, Ij C n and \I :j \ > 2, 
where for every pair Ii , Ij £ F 

Ii fl Ij = 0 or F C Ij or Ij C Ij. 

The set of all forests of n indices together with the empty forest {} is indicated by Sfj. 

For every subset / C n we indicate by Ri the operator which extracts the principal part with 
respect to aj of a multivariate meromorphic function /({o^ j-^j), where for every i,j £ I, ct,, = a/, 
and multiplies it with —1: 

Rif ■= -PP Q .lim (3.8) 

We complement this definition by setting !?{} to be the identity. 

Given all these data, we define the renormalised time-ordered product in the MS-scheme as in e.g. 
[DFKR14, Theorem 3.1] by 

Tn = (T n ) ms ■.= hm TO o ( Y, lH° T "“ } ’ (3 ' 9) 

\FG3plGF I 
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where, in the product over I £ F, Ri appears before Rj if I C J. Furthermore, for each graph F, the 
limit a = —> 0 is computed by setting a.y = ay for every i < j before taking the sum over 

the forests and finally considering the limit ar to 0. In this context we recall that, for every element 
of the sum over part of the limit ay, —> ar is already taken by applying Rj , see (3.8). 

Given the renormalised T n in the MS-scheme, the corresponding local S'-matrix may be constructed 
as 

OO • yi 

n =0 

for any local interaction Lagrangean V. 

In order to implement the minimal subtraction scheme as outlined above we first need to specify an 
analytic regularisation Hp >> of the Feynman propagator Hp on generic curved spacetimes. Afterwards 
we have to demonstrate that for all graphs T £ Q n the analytically regularised integral kernels 


appearing in 


» 


n 


Hp' J ( Xi,Xj ) 


e=(*,j)e r 


A {Hp{xi,x 3 )) h F 

l<i<j<n 


rp(a) 

n 



j2|B(r)| 


Iliey(r) ELd j S<pi(xi ) 


(3.10) 


(3.11) 


satisfy the properties necessary for the implementation of the MS-scheme. In particular we need to 
demonstrate that the distribution r^ a \ which is a priori defined only on M n \ D n , can be uniquely 
extended to the full M n without renormalisation, where the uniqueness of this extension is important 
in order to obtain a definite renormalisation scheme. Moreover, we need to show that this distribution 
Tp°‘' > £ T>'(M n ) is weakly meromorphic in a in a neighbourhood of 0, where in view of the forest 
formula it is only necessary to show that, setting = a/ for all i,j £ I, rf. is weakly meromorphic 
in ai. Additionally, we need to prove that, if it prior to regularisation is well-defined outside of 
the partial diagonal d/, then the pole of rwith ay, = a/ for all i,j £ I in ay is supported on 
di and thus local. Finally, we need to prove that our MS-scheme satisfies all properties given in 
[HW02, HW04] which a physically meaningful renormalisation scheme on curved spacetimes should 
satisfy, and we need to provide means to explicitly compute the minimal subtraction, which after all is 
the main motivation for this work. 

Our plan to construct the mentioned quantities and to prove their required properties is as follows. 


a) In Section 3.2 we construct an analytic regularisation H^ of the Feynman propagator based 
on the observation that locally Hp is of the form (2.2) up to considering instead of cr+ the half 
squared geodesic with the Feynman e-prescription ctf := a + ie. Motivated by the fact that the 
singular structure of Hp originates from the form in which crp appears, we set locally 


:= lim 


1 


e-X)+ 87T 2 \Af 2a CTp +Q 


-II- 

a 


M 2 ° 


(3.12) 


where we use the (arbitrary but fixed) mass scale M present in (2.2) also for preserving the mass 
dimension of Hp in the regularisation. 


b) In Proposition 3.7 we then prove that the relevant distributions 


i r 


n yjkyfi AM~\D n ) 

l<i<j<n u F 


(3.13) 


are multivariate analytic functions. The distribution (3.13) only displays the most singular 
contribution of r ^ (3.10), but the subleading contributions are clearly of the same form up to 
replacing some of the factors (1 + ayy) in the exponents by a^- or 0. 
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c) In order to show that t can be uniquely extended from M. n \D n to A4" in a weakly meromorphic 
fashion, i.e. that the singularities relevant for the forest formula are poles of finite order, we follow 
a strategy similar to the one used in [HW02] and consider a scaling expansion with respect to a 
suitable scaling transformation. We first argue in Proposition 3.8 that an analytically regularised 
distribution G V' (A4 n \ d n ), which can be written as a sum of homogeneous terms with 
respect to this scaling transformation plus a sufficiently regular remainder, can be extended to 
M n in a weakly meromorphic way, were the uniqueness of the extension follows from its weak 
meromorphicity. In Proposition 3.10, we give a sufficient condition for the existence of such 

(a.) 

a homogeneous expansion and we demonstrate in Proposition 3.12 that the distributions fp 
satisfy this condition. 

d) The above-mentioned results are proved by means of generalised Euler operators (see [Dal3] 
for a related concept) which can be written abstractly in terms of a scaling transformation, but 
also in terms of covariant differential operators whose explicit form can be straightforwardly 
computed as we argue in Section 3.3.1. In Proposition 3.9 we use these operators in order to 
demonstrate how the full relevant pole structure of can be computed, thus showing the 
practical feasibility of the MS-scheme. We find that our renormalisation scheme corresponds in 
fact to a particular form of differential renormalisation and expand on this by computing a few 
examples in Section 3.5. 

e) Finally, in Proposition 3.14 we prove that the MS-scheme satisfies the axioms of [HW02, HW04] 
for time-ordered products and in addition preserves invariance under any spacetime isometries 
present. 

Remark 3.1. The local Hadamard expansion (2.2) of H F and correspondingly the analytically 
continued defined in (3.12) are only meaningful on normal neighbourhoods A f of (A4, g). In order 
to define H^ and the induced distributions Tp°^(3.10) globally, we may employ suitable partitions of 
unity. Rather than providing general and cumbersome formulas, we prefer to illustrate the idea at the 
example of the triangular graph 

r r = Hpi 3 H F 23Hp 12 . := H F {xi,x 3 )H F [x 2 ,x 3 )H F {x 3 ,x 2 ) 2 

the renormalisation of which is discussed in detail in Section 3.5.3. 

We recall that a corollary of Lemma 10 of Chapter 5 in [0’N83] guarantees that there exists a 
covering C of M consisting of open geodesically convex sets such that A/) fl A/} is geodesically convex 
for every A/i, A/j G C 2 . With this C at our disposal, we define the sets 

A fi 2 := (J M x J\f C A4 2 , J\f 123 := IjA/'xA/'xA/'c M 3 . 

AAec M&c 

We call sets of the form A /12 and A /123 a normal neighbourhood of the total diagonal. This 
definition is essentially motivated by the fact that for every x G A4 we can find a normal neighborhood 
A4 G C of x in M. The squared geodesic distance a is then well defined on A/ 12 , whereas the same is 
in general not true if we replace C in the previous formula with a covering of A4 formed by all open 
geodesically convex sets. 

Setting a,j := a(xi,Xj ), we observe that < 7 12 is well-defined on A/12, and that ay 2 , 013 and <r 2 3 
are well-defined on A/123- We now consider smooth and compactly supported functions X12 G 22 (A/"i 2 ), 
X123 G 22 (A/123) which are such that Xi2 = 1 on d 2 C A/12 and X123 = 1 on d 3 C A/123. Note that 
by construction X12 and X123 vanish outside of A/12 and A/123 respectively. We may now define the 
analytically regularised distribution by setting 

4 a) := HjggHjgg X12X123 + H F , 13 H F , 2 3 H 2 F ' 12 (1 - X12) 

+ H F 13 H F23 Xi2(l — X123) > 

2 We would like to thank Valter Moretti for pointing this result out to us. 
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where the Feynman propagators are regularised as in (3.12). By construction, Tp 0 - 1 is globally well- 
defined and the analysis outlined above and performed in the following sections implies that it can be 
uniquely extended to a weakly meromorphic distribution on the full A4 3 . Moreover, the local pole 
contributions corresponding to a 12 = a/ with I = { 1,2} and 012 = 013 = 023 = Qj with J = {1,2, 3} 
are clearly independent of the choice of X 12 , X 123 and A/ 12 , A /123 such that the MS-regularised amplitude 
(rr)ms is both globally well-defined and independent of the quantities entering the global definition of 
the analytic regularisation. 

Keeping this approach to define global analytically regularised quantities in mind, we shall for 
simplicity work only with local quantities in the following. 


3.2. Analytic regularisation of the Feynman propagator Hp on curved spacetimes 

Following the plan outlined in Section 3.1, we would like to define an analytic regularisation H^' 1 of 
Hp by (3.12). To this end, we start our analysis by constructing the distribution l/cr^ +Q in A4 2 for 
a £ C \ N. As anticipated in Section 3.1 we shall make use of scaling properties of l/cr^ +a and the 
induced quantities tp°^ (3.13) with respect to a particular geometric scaling transformation. 

For every pair of points 27,07 in a normal neighbourhood M C (Ad, <7) there exists a unique 
geodesic 7 connecting x\ and Xi. We shall assume that 7 : A i-» 07 (A) is affinely parametrised and 
that 07 ( 0 ) = x\ whereas 07 ( 1 ) = 07 . For all A > 0 and all / £ V(Af n ) with J\f n C Ad™ a normal 
neighbourhood of the total diagonal d n (cf. Remark 3.1), the geometric scaling transformation we 
shall consider is 

fx ■■= A 4 ( ™ - 1 //(aq, £ 2 (A),..., * n (A)) f[ , (3.14) 

i=2 vsr) 

where g(x) is the absolute value of the determinant of the metric expressed in normal coordinates. For 
A > 1 it may happen that 07 (A) lies outside of A f n and is thus not well-defined in general. In this case 
we set f\ = 0 which is well-defined because / = 0 outside of J\f n . For later purposes, we recall that 
the determinant of the metric computed in normal coordinates centred at x\ is such that 


Vdixi) 


1 

U 2 (Xi,Xi) ’ 


where u is the Hadamard coefficient in (2.2) and u 2 is the van Vleck-Morette determinant, see e.g. 
[PPV11, (8.5)]. 

By means of this transformation, relevant information about the behaviour of a distribution in the 
neighbourhood of the total diagonal d n can be obtained. 

We recall two definitions and a few results which we shall use in the following. To this end we 
consider a distribution t in 2?'(A4” \ d n ). The restrictions of t to the sets 


O x = {(x,y) £ Ad™ | ye AT" 1 } 


for x G M are well defined if WF(f) is disjoint from the conormal bundle N*O x = {( y,k ) £ 
T*M | (fc,£) = 0,V£ £ T y O x }, see e.g. Theorem 8.2.4 in the book of Hormander [H690]. By 
assumption, these restrictions are defined everywhere on O x up to the single points O x (~1 d n . The 
extension of the restricted distributions to that point can be constructed by means of particular scaling 
properties of that distribution, see e.g. Theorem 5.2 and Theorem 5.3 in [BFOO]. However, in order 
to discuss the extension of t to the full diagonal d n it is useful to introduce the scaling degree of t 
towards d n in the following way, which is an adaption of the notion of a transversal scaling degree 
presented in [BFOO, Section 6 ] to our setting. 

Definition 3.2. Let T C T* A4™ be a closed conical set which is equal to the conormal bundle N*d n 
on d n and disjoint from N*O x for every x, and let t be a distribution in V r (X), X £ {A4™, A4” \ d n } 
i.e. t £ T>'(X) and WF(<) C T. The scaling degree of t towards d n is defined as 


sd drl {t) := inf 


£ K 


lim X w t\ = 0 

A—> 0 + 
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where the distributions t\(f) := are defined in terms of the scaling transformation (3.14), and 

where the limits are taken in the topology of V' r (A4 n ). 

If a distribution has scaling degree towards d n lower than the total dimension of the scaled 
coordinates 4 (n — 1), then it possesses a unique extension towards d n with the same scaling degree, see 
e.g. Theorem 6.9 of [BFOO]. The proof of existence of the extensions presented there is constructive 
and can be obtained as the limit to large n of t(l — 6 n ), where 9 n are suitable smoothed characteristic 
functions of a neighborhood of the diagonal whose supports become smaller and smaller for large 
n. In particular, the previous notion of scaling degree permits to control the wave front set of every 
t( 1 — 9 n ) and of the limiting distribution as well. In fact, the extended distribution can be restricted 
to every O x as well because T and N*O x are disjoint sets. The uniqueness of the extension descends 
from the uniqueness of the extension of the distributions restricted to O x and from the observation 
that the scaling degree of the restricted distributions towards the point (x,..., x) £ O x is such that 
sd(f| 0 ) < sd d n {t)- If not strictly necessary we shall drop the suffix d n in the following. Having 
this result at disposal implies that, in order to construct an extension towards d n , it is sufficient to 
construct extensions of the distribution restricted to O x . 

Remark 3.3. The distributions relevent for our analysis are the ones of the form (3.10), respectively 
(3.13). As we discuss briefly in the proof of Proposition 3.7, these distributions satisfy the micro 
local spectrum condition introduced in [BFK96] and, consequently, also the wave front set condition 
required in Definition 3.2. 

The scaling degree towards a partial diagonal may be defined in analogy to Definition 3.2. The 
same geometric transformation (3.14) can be used to introduce relevant homogeneity properties of a 
distribution. 


Definition 3.4. A distribution t £ V\M n ') or t £ T>'(A4 n \ d n ), which satisfies the equality 

= VA > 0 


under transformations of the form (3.14) for all / £ T)(N n \d n ) and for a S £ C, is called homogeneous 
of degree 5. 


These definitions imply by direct inspection that a distribution which is homogeneous of degree 
S and whose wave front set restricted to d n is contained in N*d n has scaling degree —Re(<5). We 
further notice that homogeneous distributions t £ T>{A4 n \ d n ) whose WF(t) fl A4 is contained in 
N*d n , possess unique extensions to M. n with the same degree of homogeneity S if — (5 + 4(n — 1)) ^ N. 
A proof of this claim can be obtained following very closely the proof of Theorem 3.2.3 in [H690]. 
In particular, the uniqueness of these extensions descends once again from the uniqueness of the 
extensions of the distributions restricted to A4 X , whereas the existence of the extension is obtained 
employing a fundamental solution of the Euler orperator E\ we shall introduce below in (3.17). 

In the next proposition we introduce the distributions we shall use as building blocks for the 
construction of regularised Feynman propagators on Lorentzian manifolds. Although not completely 
analogous, the construction we are going to present is similar to the extension of Riesz distributions to 
curved spaces presented in Section 1.4 of [BGP07]. In particular, here we shall discuss the boundary 
value of l/(cr + ie) a for e — > 0 while ordinary Riesz distributions are related to the antisymmetric part 
of a different boundary value of the functions l/a a . 

Proposition 3.5. Consider a normal neighbourhood A4 C A4 2 of c ?2 (cf. Remark 3.1) and the 
following expression for a £ C and f £ VijKff) 




:f(x, y)dp g (x)dp g (y). 


er£’ /' e->o+ (a(x,y) + ie)' 

Then the following statements hold. 

a) l/<Jp restricted to V{J\f 2 \ cfe) is a distribution which is weakly analytic in a. 
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b) l/<Jp is homogeneous of degree —2a with respect to transformations of the form (3.14) for 

f £ T>{M 2 \ d 2 ). 

c) l/crp is well-defined as a distribution on M 2 for 2a — 4 ^ N. Furthermore, for all f £ 2?(A/ 2 ) 
(1/cTp, f) is analytic for 2a — 4 ^ N and meromorphic for a £ C with simple poles at 2a — 4 £ N. 

Proof, a) For every x £ Mi we fix a normal coordinate system £ x : y — > M 4 in order to parametrise 
points y in a normal neighbourhood of x. Consequently, on M 2 the squared geodesic distance divided 
by 2 can be easily expressed as 

<?{x,y) = ^v(€x(y),€x(y)) = 

where 77 is the standard Minkowski metric given in Cartesian coordinates. Furthermore, 




'M JR 4 


(£S£za + * e )' 


:f(x,€x)Vg(€x) d A f, x dfigix). 


(3.15) 


which is well defined for f £ V(J\f 2 ). 

Observe that l/(£ a £a)“ for £ a £ {z £ C 4 | Im (z) £ V ± }, where V ± is the forward or past light 
cone with respect to the Minkowski metric, is analytic both in £ and a. Furthermore, in the limit 
e -+ 0+, l/(£°£a + ie) a can be seen as the boundary value of that analytic function. Since this function 
grows at most polynomially for large 1/Im(£ a £ a ) its boundary value defines a distribution, see e.g. 
[H690, Theorem 3.1.15]. The analytic dependence on a is weakly preserved in the limit e —> 0 + , and 
thus the resulting distribution is weakly analytic. 

b) The transformation defined in (3.14) acts on points parametrised by normal coordinates as £ —> A£. 
Furthermore, l/(£ a £o)“ on A C C 4 is homogenous of degree 2a with respect to the transformation 
£ —»• A£. The statement follows from this observation, taking into account (3.14) and (3.15). 

c) Theorem 3.2.3 in [H690] ensures that the distribution t = l/(£ a £ a )“ £ ^'(R 4 \ 0) has a unique 
extension t to 0 preserving the degree of homogeneity for every 2a — 4 ^ N. Hence, I<£> t defines 
a distribution on M x R 4 . Finally, notice that there exists a neighbourhood O of the diagonal in 
Ad 2 where f(x,£ x )y/g(£ x ) is a smooth and compactly supported function for every / £ Cfi°((D). 
Consequently, the statement follows by choosing a partition of unity adapted to O. We refer to the 
discussion after Definition 3.4 for general arguments why in fact Theorem 3.2.3 in [H690] can be 
generalised to extensions of suitable homogeneous distributions to the diagonal rather than a point. □ 


The previous proposition guarantees that 1 /<jp is weakly meromorphic in a with simple poles at 
2a — 4 £ N. This property is preserved under taking linear combinations and multiplication by smooth 
functions. Consequently, the analytically regularised Feynman propagator defined by (3.12) is 
well-defined on a normal neighbourhood of the diagonal and weakly meromorphic in a. 

Proposition 3.6. Consider a normal neighbourhood M 2 of the diagonal d 2 £ Mi 2 . The following 
statements hold for the analytically continued Feynman propagator H^ £ V(M 2 ) defined in (3.12). 

a) lim = II F . 

b) WF{H ( p ] ) C WF(£T f ). 

c) The scaling degree of tends to —00 when the real part of a tends to 00 . 

Proof. The proof of this proposition follows from the properties of ct^ +q obtained in Proposition 3.5. 
In particular, a) and c) can be directly obtained from the weak analyticity, while b) follows from 
the fact that the distribution l/(Jp is well defined on M 2 \ <^2 where it coincides either with 1 /<r“ or 
with 1 /cr“. In order to analyse the wave front sets of 1 /a± , we pass to a normal coordinate system 
and obtain 1 /<t“ = 2/(£ a £ a ± ?'e£°) a . This distribution can be extended to a tempered distribution 
for every a and thus its Fourier transform can be directly computed. One finds that for 1 /cr±, 
only the null future/past directed directions do not decay rapidly, consequently Fd^ restricted to 
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J\f 2 \ d 2 has WF(U^) C WF (H F ). Finally, we observe that the extension of to J\f 2 may possess 
further singularities supported on the diagonal with singular directions orthogonal to d 2 . Hence, 
WF (H ( F a) ) c WF (H f ) still holds for H < a) £ V'(U 2 ). □ 


We are now able to discuss the analytical regularisation r^ (3.10) of the distributions Tr given in 
(3.6) which appear in the graph expansion (3.4) of the time-ordered products T n (3.1). As anticipated 
in Section 3.1, owing to the form of H given in (3.12) the relevant distributions which need to be 
discussed are introduced in (3.13) and analysed in the following proposition. 


Proposition 3.7. The operation 


(tr a) j) ~ [ II —7- (i+ a ) f dxi ■ ■ ■ dx * 


(3.16) 


defined for f £ V(A4 n \ D n CiAf) where J\f is a normal neighbourhood of the total diagonal (cf. Remark 
3.1) has the following properties. 

a) is distribution on A4 n \ D n D A/". 

b) (t^, is a continuous function for a = {aij}i<j £ M n ( n-1 )/ 2 . 

c) It p“\ is analytic for every ctij with i < j and thus a multivariate analytic function. 


Proof, a) The domain A4 n \ D n n TV" is a disjoint union of connected components. On every connected 
component C a F (xi, Xj) equals either cr+(xi, Xj) or <j+(xj, xf) depending on the causal relation between 
Xi and Xj which is fixed in C. Hence, on C, the wave front set of cr F (xi, sq) -1 is contained either 
in V+ or V_, where V + /_ = {(x,x',k,k') £ T*M 2 \0 ,(x,k) ~ (x', — k'), k < / >0}. Consequently, 
a F (xi,Xj ) _1 satisfies the Hadamard condition up to a permutation of the arguments. The very same 
holds for the distributions cr F (xi, for every hj and every which have been discussed in 

Proposition 3.5. 

Owing to the form of their wave front set, the pointwise products of these distributions present in 
<p“^ are well-defined because the Hormander-criterion for multiplication of distributions is satisfied. In 
fact, up to some fixed permutation of the arguments (aq,... x n ), t^ satisfies the micro local spectrum 
condition introduced in [BFK96]. Hence t^ is a well-defined distribution on every connected 
component C of M n \ D n n A f and thus it is well-defined also on M n \ D n D Af. 

b) In order to check continuity for a = {a?,j}i<j £ R n (" -1 )/ 2 in a fixed point a we may analyse the 

(cx) 

distribution on a fixed connected component C of the domain of t r and factorize the distribution in 


two parts. In fact, due to the wave front set of fp“' ) on C the factorisation £p°^ = t\f* 


id) ■ 

Tp is unique 


-. For /3 in a sufficiently small neighbourhood 


where the integral kernel of r i 13 ' 1 is — 7 —— 

1 a F (Xi ,Xj, - 

of 0, t^P' 1 is an integrable function which is differentiable for /3 = 0 as can be obtained by dominated 
convergence. Finally, the continuity is preserved by pointwise multiplication with fp Q \ 

c) For an arbitrary but fixed pair of indices i,j, aij appears in the product displayed in (3.16) as 
l/cr F (xi, Xj) li:i ( 1+ai ^ and we have already analysed the analyticity property of such a distribution in 
Proposition 3.5. We shall thus interpret t as a composition of distributions, namely as l/crpf 3 o z 
where z is an operator which maps V(M n \ D n f)N) to V'(M 2 \ D 2 nJ\f 2 ) for a suitable Af 2 D D 2 = d 2 . 
The e-regularised integral kernel of z corresponds to the product present in (3.16) with the factor 
removed. Because of the singular structure of z, for every / € V(M n \ D n fI A/"), (z, /) is in 
fact a compactly supported smooth function supported on M 2 \ D 2 fl J\f 2 . Hence, the analysis of its 
composition with l/cr^.’ 3 is straightforward. These considerations imply separate analyticity of t < p*' 1 in 
each whereas joint analyticity follows from the continuity proved in b). □ 
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3.3. Generalised Euler operators and principal parts of homogeneous expansions 

The next step in the strategy outlined at the end of Section 3.1 is to extend the distributions which 
are a priori defined only outside of the union of all partial diagonals D n in a normal neighborhood N of 
the total diagonal (cf. Remark 3.1) to D n n N and to show that this extension is weakly meromorphic 
in ai upon setting = oq for all i,j £ / C {1,... ,n}. As anticipated, we shall prove this by 
using particular homogeneity properties of with respect to the scaling transformations (3.14). 

(a.) 

Even if tf is not homogeneous in the strong sense of Definition 3.4, it has weaker homogeneity 
properties which are still strong enough in order to obtain the wanted results. In this section we analyse 
analytically regularised distributions satisfying this weaker homogeneity condition, provide sufficient 
conditions for this weaker homogeneity to hold and show how the principal part of a distribution of 
this type can be efficiently computed. 

To this avail, we consider a normal neighbourhood N n of the total diagonal d n (cf. Remark 3.1) 
and define the generalised Euler operator E p : V{N n ) —> V{N n ) by 


E P f(xx n ) := (-l) p (a- 4 ^- 1 )/,(*)) 


(3.17) 


where the scaling transformation (3.14) is used. We then consider a family of distributions t £ 
V' r (Af n \ d n ) C V{N n \ d n ) (cf. Definition 3.2) defined for a in some neighbourhood O of 0 £ C and 
assume that t^ can be expanded as 


t^ = J2 t{ k a) +r (a) - 

k—0 


where are homogeneous with degree a*, = — S a + k whose real part is smaller or equal to —4(n — 1) 
and a remainder £ V {N n \ d n ) which has scaling degree smaller than 4 {n — 1) and can thus 
be uniquely extended to d n for every a £ O by [BFOO, Theorem 6.9], see also the discussion after 
Definition 3.2. Owing to its homogeneity, every t^ can be rewritten by means of the generalised 
Euler operator E p as 



1 

]lj=d(«fc+j + 4(n 



(3.18) 


Note that, E p f{x i,...x n ) is smooth and vanishes for y = (x±,.. .x n ) —i x = (aq,..., aq) as 
C\y — x\ p , i.e. it is in the class 0{\y — a;| p ). For this reason, if p is chosen sufficiently large as 
p > —ak — 4(n — 1), oE p possesses a unique extension to d n . We recall that, in order to renormalise 
t(“) for a = 0 in the MS-scheme, we have to subtract its principal part before computing the limit of 
vanishing a 

((4)ms,/) := lim ((4“ } >/) ~Pp(4 a) >/)) ■ 

However, if we use the representation of provided by the right hand side of equation (3.18), its poles 
are manifestly exposed and can be easily subtracted. We recall that, since the original distribution 
is well defined on A f n \ d n even for a = 0, the principal part we are subtracting can only be supported 
on d n . We summarise this discussion in the following proposition. 


Proposition 3.8. Consider a normal neighbourhood Af n of the total diagonal d n and a distribution 
t £ V[J\f n \ d n ). Assume that t^ £ V{N n \ d n ) is an analytic regularisation oft, i.e. is weakly 
analytic for a in a neighbourhood O of the origin of C and lim^-m t^ = t. Moreover, assume that 
can be decomposed as 

m 

t {a) = J2 4 a) + r(a) 

k =0 


where t^ are weakly analytic distributions which scale homogeneously under transformations of the 
form (3.14) with degree ak = —<+ k and r^ is a weakly analytic distribution whose scaling degree 


18 




towards d n is strictly smaller than 4(n — 1). Then, ifWF(t) and WF(ff“)) do not intersect the 
conormal bundle of the orbits of the scaling transformation (3.14), the following statements hold. 

a) tcan be extended to £ V'(Af n ) for every a £ O \ {0}. 

b) *(“) is weakly meromorphic for a £ O with possible poles for a = 0 and it is the unique weakly 
meromorphic extension of t^. 

c) The pole of t^ in 0 is supported on d n . 

d) The limit a —> 0 can be considered after subtracting the pole part, namely 

(tms,f) ■= lim ((i (Q) ,/) -PP 

is well-defined for all f £ T>(N n ) and t ms is an extension oft which preserves the scaling degree. 

Proof. The proof of a) and b) follows from the discussion after Definition 3.4 as an application of 
[H690, Theorem 3.2.3] to every restricted to every orbit of the action of the scaling transformation 
(3.14). Furthermore, since the scaling degree of r ^ towards d n is strictly smaller than 4 (n — 1), r 
possesses an unique extension towards d n , cf. [BF00, Theorem 6.9]. 

In order to prove c) we note that the original distribution t^ defined on Af n \ d n is weakly analytic 
and that an explicit construction of the weakly meromorphic extension to A f n is provided by (3.18), 

choosing for every component t a sufficiently large p and using again [BF00, Theorem 6.9]. Hence, 
the poles of t^ can only be supported on d n . For this reason, after subtracting the principal part of the 
distribution the limit a —> 0 can be safely taken. The such obtained distribution prior to considering 
the limit a —> 0 coincides with t^ on Af n \ d n and the same holds in the limit a —> 0. Consequently 
f ms is an extension of t. Finally, sd(f ms ) = sd(t), because our assumptions and the above analysis 
imply that sd(i^) = sd(^ a ^) = Re (<5o,), sd(pp(^“^)) < Re ( S a ) and lim^o Re (<5 Q ) = sd(f). □ 

We now discuss how equation (3.18) can be used in order to regularise the most singular part of a 
distribution t^ which is known to be of the form t^ + r ^ but where the distributions 

t^ are not explicitly known. To this end, observe that equation (3.18) implies 

m fp —1 \ 

(fW,F; p /) = ^ Y[(a k +j + 4(n-l)) (t^\f) + (r^,E p f). 

k=0 \j=0 J 

Moreover, we may assume without loss of generality as in Proposition 3.8 that the homogeneity degrees 
a k of t^ 1 are of the form a k = —S a + k where Re ((5 a ) is the scaling degree of t^ a \ Consequently, 
tg 0 ) is the contribution with the highest scaling degree which may be extracted by introducing the 
coefficients 

p -1 

Cfc := (a k +j + 4 (n - 1)) 
j=0 

and considering 


(t {a \E p f) - co /) = - c o) (4 Q) »/) + (r (a \E P f) - co (fW, /) , (3.19) 

k =1 

where the distribution on the right hand side has a scaling degree smaller than Re (c> Q ) = —Re (ao). 
Hence, although in general the distribution t^ does not scale homogeneously, equation (3.18) still 
holds up to distributions with a lower scaling degree. Knowing the decreasing degree of homogeneity 
of the components in the expansion of t^ a \ we may use a recursive procedure in order to expose the 
pole part of this distribution. In fact, the previous discussion straightforwardly implies the validity of 
the following proposition. 
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Proposition 3.9. We consider a distribution with the properties assumed in Proposition 3.8 and 
set 

u 0 ■■= u k +i ■■= c k u k - Uk ° E Pk , 0 < l < m 

where p k are the smallest natural numbers chosen in such a way that p k + Re (a k ) + 4(n — 1) > 0 
and c k := YVU (afc + j + 4(n — 1)). Then, in order to expose the poles of t^ a \ we may invert the 
recursive definition of u k obtaining 

t = — ( uo° E po H- (u\ o E P1 + ■ • • H-( u n o E Pri + u n +i) J J . (3.20) 

Co \ ci \ c n j j 

In order to be able use the previous results for our purposes, we provide in the next proposition a 
criterion which is sufficient to ensure that a generic distribution can be decomposed into the sum of a 
homogeneous distribution and a remainder with lower scaling degree. We shall use this criterion in 

(ck) 

order to prove that the distributions t r defined in Proposition 3.7 have the desired property. 

Proposition 3.10. Let A f n be a normal neighbourhood of the total diagonal d n and suppose that 
t € V\N n ) has scaling degree si towards d n under transformations of the form (3.14) and that there 
exists an a with —Re (a) = si such that t o + a + 4(n — 1)) has scaling degree S 2 < s\. Then t can 
be decomposed into the sum of a homogeneous distribution with degree a and a remainder with scaling 
degree smaller than or equal to S 2 - 

Proof. We start by observing that, for every test function / G T>(Af n ), 

F(XJ) := (t,f 1/x ) 

is a continuous linear functional of / which is smooth in A for A > 0. Moreover, since the scaling 
degree of t is si, X a F(X, f) vanishes in the limit A —> 0 for every a > Si and for every / G T>(M n ). Let 
us now consider 

G( A, /) := ((—.Ei + a + 4 (n — l))t, /i/a) ■ 

G(A, •) is again a family of distributions on J\T n which depends smoothly on A for positive A. Furthermore, 
A“G(A,/) vanishes in the limit A —>■ 0 for every a > s 2 and every / G T>(N n ). Hence, A“ _1 G(A, •) 
tends to 0 in V(No) for A —> 0 and, additionally, the Banach-Steinhaus theorem implies that 

|A a G(A,/)|<G^|9“/|, (3.21) 

ot<k 

for every a > S2, uniformly for / supported in a compact set K C A f n and for suitable C and k which 
do not depend on A. 

After these preparatory considerations, we observe that G and F are related by the generalised 
Euler operator in the following way 

G(A,/) = A-“ +1 ^A“F(A,/). 

We can invert this relation to obtain 

F (A, /) = ^ j Q X ~X a ~ 1 G(X, f)d~x , 

where, C(f) is a suitable constant which depends on /. We want to prove that G(-) is in fact a 
distribution. To this end, we note that, owing to the bound (3.21), the integral in A can be performed 
and the result of this integration is a distribution for every A > 0 because Re (a) > S 2 - This implies 
that 

C(f) = F(l, f) — f X a ~ 1 G(X, f)dX 
Jo 

is a distribution because it is a linear combination of distributions. By construction E(l, /i/a) = F( A, /) 
and G o [E\ + a) =0, hence G is a homogeneous distribution of degree a. By means of a direct 
computation we also find that the scaling degree of the remainder F( 1, /) — G(/) is smaller than or 
equal to the scaling degree of G which is S 2 - □ 
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3.3.1. The differential form of generalised Euler operators and homogeneous expansions 
of Feynman amplitudes 

In order to make the previous discussion operative, we have to analyse the action of the generalised 
Euler operators E p appearing in (3.18) on test functions. In fact, we shall see that E p corresponds to a 
particular geometric partial differential operator. To this end, we observe that E p = ( E\ — (p— l))E p _i. 
Hence, knowing the differential form of the generalised Euler operator E\, it is possible to construct 
recursively every E p . 

Regarding the differential form of E\, we note that it can be written in terms of the geodesic 
distance and the van Vleck-Morette determinant 3 * u 2 as 

n 

Eif(x u ...,x n )=^2 (a a (xj)V^ - (2a a ( Xj )V^ log(u(xj, an)))) f(x i,..., x n ), 

3 = 2 

where Va 3 indicates the a-th component of the covariant derivative computed in Xj and cr a (xj) := 
V X: > a cr(xi,Xj). Considering the adjoint E£ of E pi we have t o E p = E^t where, using the relation 
□cr + 2er a V a log(zt) = 4, we find for p = 1 

n 

E\t{x 1 , ...,X n ) = ( ~ V a 3 ^(Xj) ~ ^(Xj) (Va 3 log(u(Xj, Zl)))) t(x 1 , . . . , X n ) 

1=2 


= - ( 4(n — 1) + ^ a a (xj)V Xj ] t(x u ..., x n ). 

i=2 


(3.22) 


We finally observe that the recursive identity for E p implies that also E | can be constructed recursively 
starting from E[ as 

= £!_,(£} _ (p _ i)). (3.23) 

We proceed by showing that upon applying E\ introduced in (3.22) to a distribution t p Q) of the 
form 


L r ~ 


n 


„ ^ ap(xi, Xj) li j( 1+OLi i) 

which has scaling degree sd(fp°^) = Yhi<j 2Zij(l + Re (a*j)) towards the thin diagonal d n , the result is a 
term proportional to plus a remainder which has lower scaling degree as foreseen in (3.19). Hence, 
Proposition 3.10 implies that can be written as a homogeneous distribution plus a remainder with 
lower scaling degree. If the scaling degree of the remainder is not sufficiently low, we reiterate the 
procedure in order to obtain a full almost homogeneous expansion of the desired form. 

In order to analyse this issue we shall only consider the relevant differential operator on A4 n 
appearing in E[, namely, 


P ; = - 

1=2 


(3.24) 


We start by analysing the action of p on cr(£ 2 , £ 3 ) f° r x^,x^ in a normal neighbourhood of the point 

Xl- 

Lemma 3.11. Let J\f Xl be a normal neighbourhood of the point X\ and let X2,X3 £ Af xi - Then, 

pcr(x 2 , x 3 ) = 2a(x 2 , x 3 ) + G(xi,x 2 , x 3 ) 

where G is a smooth function which vanishes in the limit x 2 , x 3 —> X\ as a monomial of order 4 in the 
normal coordinates of x 2 and x 3 centred in X\. 


3 Recall that the square—root of the van Vleck—Morette determinant coincides with the Hadamard coefficient u 

appearing in (2.2). 
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Proof. Using the notation in the proof of Proposition 3.5 we write the action of p on 023 := cr{pc 2 ,^ 3 ) 
as 

P °23 = fa{X2)v%3 + £b> (% 3)^3 . 

Recall that crf 3 is the covector in T* o M cotangent to the unique geodesic joining X 2 and ,t 3 , that — cr\ 3 
is equal to the parallel transport of crf 3 from X 2 to x 3 along the geodesic 7 joining the two points, and 
that £ c (a 7 ) := a c (xi,Xi). 

Let us parametrise the image of 7 with an affine parameter A such that cc(0) = X 2 and x(l) = x 3 . 
In order to simplify the notation, we indicate by t( A) the tangent vector of the geodesic in x(A). As 
argued before, we have 

t a (0) = ct% 3 , and t b '(l) = -a^. 

Consequently, 

po -23 = £“ta(0) -£ b t b (l) = -J f(£ a t a )(A)dA 

= - f t a X7 t f a dA = [ t a t b <T ab (x(X),x 1 )dX , 

Jo Jo 

where a ab := V Q Vhcr. If we now consider the covariant Taylor expansion of cr ab (x(A), xf) around cc(A) 
(see e.g. [PPV11]), we find that E ab (x,x 1 ) := a ab (x,x 1 ) — g a b(x) is a smooth function that vanishes 
for x —> X\ as 0(a(x,xi)), hence 

po 23 = [ t a t b g ab (x(A))dA + f t a t b E ab (x(X),xi)dX = 2a(x 2 ,x 3 ) + G(x 1 ,x 2 ,x 3 ), 

Jo Jo 

where the remainder is smooth because of the smoothness of the metric g and can be further expanded 
as 


G(xi,x 2 ,x 3 ) = / t a (X)t b (X) (a ab {x{X), xi) - g ab (x(X))) dX (3.25) 

Jo 

t a (X)t b (X)t c (X)t d (X)R acbc i(x(X))dX + ■■■ = 0(|£(z 2 )| 4 + |^(s 3 )| 4 ), 

where the absolute value of the normal coordinates |£(a;j)| of 27 , i = 2, 3 is intended in the Euclidean 
sense. □ 



We are now in position to analyse the action of p on the distribution introduced in (3.16). 

Proposition 3.12. The distribution t^ introduced in (3.16) can be written as a sum of homogeneous 
distributions with respect to scaling towards the total diagonal d n plus a remainder. The degrees of 
homogeneity of these homogeneous distributions are contained in the following set 


k — N ) 2lij(\ T 01 ij)i k G N U {0} 

l<£<j<n 

Proof. We perform this analysis with e in ap taken to be strictly positive. We start by applying p 
given in (3.24) to t!'"' 1 . Thanks to the results stated in Lemma 3.11 we have 

= ct ( r a) + 4 a) , 

where the constant C is 

C= 2lij(l + aij). 

l<£<j<n 
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Furthermore, Lemma 3.11 and in particular (3.25) implies that the remainder has a scaling degree 
towards d n which is lower than the one of t [P* by at least two, 

sd(rp“ } ) < sd(4 a) ) - 2 = ^ 2Z W (1 + Re (ay)) - 2 • (3-26) 

l<i<j<n 


(ot) 

Proposition 3.10 then implies that the distribution ip ; can be written as a homogeneous distribution 
of degree C plus a remainder with lower scaling degree. 

In order to finalise the proof we need to control the recursive application of p, therefore we discuss 
the application of p on p n t T for an arbitrary n. Let us start with n = 1. In this case, we observe 

(a.) 

that the relevant contribution is the one given by the remainder pr-^ , which reads 





G{xi,Xj,Xj) ( a ) 

a F (xi,Xj) r 


Note that for every i < j, a F (xi,Xj)t^ has the same structure like fp“\ but the scaling degree 
sd(tp ) + 2, whereas G(x\,Xi,Xj) defined in (3.25) is a smooth function whose Taylor expansion for 

(ck) 

Xi,Xj around Xi starts with components of order 4. Hence, if we apply p to rf. we obtain a constant 
multiple of rjP* plus a remainder which has scaling degree lower or equal to sd(?’p Q *) — 1, where the 
difference with respect to (3.26) stems from the fact that G can be expanded as a polynomial in a a (xi) 
whose lowest components are monomials of degree 4 multiplied by curvature tensors. These monomials 
are homogeneous and thus contribute to the degree of homogeneity of prf, , while the contributions 
in G with degree higher or equal to five influence the scaling degree of the remainder. Repeating 
this analysis for a generic n, we find that similar results hold when p is applied recursively to the 
remainder. 

Consequently, an iterated application of Proposition 3.10 implies that the distribution can be 
written as a finite sum of homogeneous distributions plus a remainder. Furthermore, since the scaling 
degree of these distributions is always finite, the degree of homogeneity of these components is finite 
as well. Q 

As outlined at the end of Section 3.1, we can use Proposition 3.12 in conjunction with the 
propositions 3.8 and 3.9 in order to extend the distributions ip ' in a unique and weakly meromorphic 
fashion to the union of all partial diagonal D n in a normal neighborhood of the total diagonal (cf. 
Remark 3.1) and in order to compute the relevant pole part of this extension as used in the forest 
formula, cf. (3.9), (3.10) and (3.12). To this avail, we stress that Proposition 3.12 holds in particular 
for any subgraph T/, I C {1,..., n} of T and the corresponding distribution which is obtained 

by omitting all factors in t r ' which correspond to edges not contained in T/. Finally, the recursive 
structure of the forest formula (3.9) implies that we are not dealing only with expressions of the form 
\aij=a I Vi,jeh but also with expressions which are of this form up to a subtraction of their principal 
part. However, our above analysis and in particular the discussion in the proof of Proposition 3.8 
implies that the propositions 3.12 and 3.9 also hold in this case. 

Remark 3.13. Proposition 3.9 and the above analysis imply that our renormalisation scheme is in 
fact a particular form of differential renormalisation. Notwithstanding, the advantage of formulating 
this scheme in terms of analytic regularisation and minimal subtraction is the ability to define the 
renormalisation scheme in a closed form at all orders by means of the forest formula (3.9). 

3.4. Properties of the minimal subtraction scheme 

We conclude the general analysis of the renormalisation scheme introduced in this work by demonstrating 
that this scheme satisfies - up to one property we shall mention at the end of this section - all axioms 
of [HW01, HW02, HW04] which, as argued in these works, any physically meaningful scheme to 
renormalise time-ordered products should satisfy. We refer to these works for a detailed formulation 
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and discussion of these axioms. In addition to showing these properties of the scheme, we also argue 
that it preserves invariance under any spacetime isometries present. 

Proposition 3.14. The time-ordered product T n defined by means of (3.9), where the quantities 
appearing in this formula are defined by means of (3.8), (3.11), (3.10) and (3.12), and were we recall 
Remark 3.1, have the following properties. 

a) Tn is symmetric and satisfies the causal factorisation condition. 

b) T n is unitary. 

c) T n is local and covariant. 

d) T n satisfies the microlocal spectrum condition. 

e) Tn is (f-independent. 

f) T n satisfies the Leibniz rule. 

g) T n satisfies the Principle of Perturbative Agreement for perturbations of the generalised mass 
term p in the free Klein-Gordon equation Pcf := (—□ + p)(j> = 0. 

h) If the spacetime (A i,g) has non-trivial isometries and if the Feynman propagator Hp is chosen 
such as to be invariant under these isometries, then T n is invariant under these isometries as 
well. 

Proof, a ) holds because we constructed the renormalised time-ordered product by means of the forest 
formula (3.9) and because, as implied by Proposition 3.8, all counterterms subtracted in the forest 
formula are local. 

b) Unitarity holds because the operation of extracting the relevant principal part of a regularised 
amplitude t^ commutes with complex conjugation (even if a is not real). 

c) The regularised amplitudes Tp 0 "* satisfy locality and covariance. Upon setting = a/ for 
i,j £ I C {1,..., n}, Tp Cl ' ) is weakly meromorphic in a/. Thus locality and covariance holds for each 
term in the corresponding Laurent series and consequently also after subtracting the principal part of 
this series. 

d) As argued in the proof of Proposition 3.7, the distributions fp 0 " 1 defined in (3.16) satisfy the 
microlocal spectrum condition, i.e. they have the correct wave front set. Consequently, the regularised 
amplitudes Tp“^ have the correct wave front set as well. As Tp“^ is weakly meromorphic in the sense 
recalled in the proof of c), each term in the corresponding Laurent series has a wave front set bounded 

(a) 

by the wave front set of Tp . Consequently the microlocal spectrum condition holds after subtracting 
the principal part and considering the limit of vanishing regularisation parameters. 

e) This property follows directly from the construction. In particular the subtraction of counterterms 
is defined in terms of numerical distributions and independent of the field f>. 

/) In analogy to b), the Leibniz rule holds because the operation of extracting the relevant principal 

(a.) 

part of a regularised amplitude Tp commutes with all partial differential operators. 

g) The Principle of Perturbative Agreement for perturbations of the generalised mass term p 
demands essentially that upon setting p = po + pi, the renormalisation of T n commutes with the 
operation of perturbatively expanding quantities in p\ around po. A Feynman propagator Hp depends 
on p only via the Hadamard coefficients v and w in (2.2). However, in the definition of the analytically 
regularised Hp in (3.12) and the corresponding regularised amplitudes Tp Q) defined in (3.10), these 
coefficients are not altered but only the <r-dependent terms multiplying these coefficients are modified. 
Consequently, the analytic regularisation and minimal subtraction scheme we consider commutes with 
a perturbative expansion in p\ around po- 

h) As recalled in g) all operations in our analytic regularisation and minimal subtraction scheme 

act directly on quantities defined entirely in terms of the geometric quantity er. As cr is invariant under 
any spacetime isometries present, the renormalisation scheme preserves this invariance. □ 
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Remark 3.15. Note that the Principle of Perturbative Agreement (PPA) as introduced in [HW04] 
also poses conditions on 71, i.e. the renormalisation of local and covariant Wick polynomials, which we 
omitted in our analysis, cf. Footnote 1 on page 10. However, given 71 for n > 1, 71 can be adjusted 
in order to satisfy the PPA for changes of /.i by using e.g. [DHP15, Theorem 3.3]. Moreover, the 
PPA as introduced in [HW04] further demands that, setting g = g 0 + gi, the renormalisation also 
commutes with perturbatively expanding quantities in g\ around an arbitrary but fixed background 
metric g 0 - Since a depends on g, it is not easy to check whether a perturbative expansion in gi 
commutes with our analytic regularisation and minimal subtraction scheme and thus it might well be 
that the renormalisation scheme discussed in the present work fails to satisfy this part of the PPA. 
However, if this is the case, the scheme can be modified according to the construction in [HW04] in 
order to satisfy also this condition while preserving the other properties in Proposition 3.14, including 
the invariance under any spacetime isometries present. 

Remark 3.16. We have omitted the explicit dependence of renormalised quantities on the mass 
scale M appearing in the analytically regularised Feynman propagator H^’ ] (3.12), but our analysis 
implies that the dependence of these quantities on M is such that all renormalised quantities are 
polynomials of (derivatives of) log ( M 2 a F (xi,Xj )), see also the examples in the next section. Thus, 
the renormalisation group flow with respect to changes of M may be easily computed. 

3.5. Examples 

In this section we illustrate the method developed in Section 3.3 to explicitly compute renormalised 
quantities in our scheme by considering first the example of the fish graph and the sunset graph, i.e. 

for n = 2,3. These pointwise powers of the Feynman propagator are the only ones occurring 
in renormalisable scalar field theories in four spacetime dimensions. Afterwards we will consider a 
triangular graph in Section 3.5.3 in order to illustrate the method in the case of more than two vertices. 
Recalling Remark 3.1, we shall work only on subsets of the spacetime where the geodesic distance is 
well-defined without loss of generality. 

In the special case of AJ., we are dealing with distributions which are already defined on M 2 \ d ,2 
and have to be extended to A4 2 . In order to accomplish this task we shall use (3.18) in order to 
expose the poles before subtracting them. In this context, we note that E\ given in (3.22) applied to a 
distribution t whose integral kernel t(a F ) depends on x, y only via a F (x,y), can be further simplified. 
In particular, introducing ti(crp) such that \7 a ti(cr F ) = a a t(cr F ), we have 

E\t(a F ) = - (4 + u a V a ) ^u) = -\/ a a a t-2a a (\/ a log(u))t(a F ) (3.27) 

\/ 7/ 

= -□t 1 (af)-2^V Q f 1 (a F ), 

u 

where x is considered to be arbitrary but fixed and all the covariant derivatives are taken with respect 
to y. 

3.5.1. Computation of the renormalised fish and sunset graphs in our scheme 

We recall that the Feynman propagator A F (x,y) := (4>(x) -t a <t>(y)} n in any Hadamard state f 1 is 
locally of the form 

A F {x,y)= ( u ^' v ^ +v(x,y)log(M 2 a F (x,y))) +w(x,y), a F :=a + ie. (3.28) 

8n z \a F (x, y) ) 

From (3.28) we can infer that, in order to renormalise A|, and A^,, i.e. in order to extend them from 
M 2 \d 2 to A4 2 , we need to renormalise the three distributions 

1 


because all other occurring powers of a F 
scaling degree for y —> x smaller than 4, and thus can be uniquely extended to the diagonal. To this 


log ( M 2 a F ) 


1 

73 " ’ 


(3.29) 


v F v F 

i.e. <jp m log"(up) for m € {0,1} and n £ {0,1, 2, 3} have a 
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avail, we define 


a a !-a n := V a „ • ■ • V ai a [B]{x) := B(x,x ), 

where the covariant derivatives are taken with respect to x and B is a general bitensor, and recall the 
following basic identities satisfied by er: 

(T a (T a = 2 a , cr a bcr b = cr a , Da = 4- 2 ° ^ aU . (3.30) 

u 

For our purposes, it will prove useful to use the last identity in the form 

Oap=4: + faF with f :=—2— -— , (3.31) 

ucjf 

where / is a distribution, which, considered as a distribution in y for fixed x, has scaling degree zero for 
y —> x as can be seen from the covariant Taylor expansion u = [it] + ([V 0 it] — V a [it]) er° + 1Z U = 1 + 1Z U , 
where the remainder 1Z U vanishes towards the diagonal faster than a a (see e.g. [PPV11, Section 5]). 

Remark 3.17. As / has vanishing scaling degree for y —> x, the pointwise product f(x,y)t(x, y) with 
any bidistribution t of scaling degree for y —> x lower than 4 may be uniquely extended to the diagonal. 
However, we will also encounter expressions which are naively of the form f(x,y)5(x,y) and which 
are a priory ill-defined because / is in general divergent for x and y light-like related, and thus not 
continuous on the diagonal. Notwithstanding, the distribution f(x,y)5(x,y), which is well-defined 
and identically vanishing outside of the diagonal x = y, may be extended to the diagonal. In fact, our 
scheme, in which expressions of the form f(x,y)8(x,y) appear as a —> 0 limits of particular weakly 
analytic expressions, provides a unique and non-vanishing extension of /( x, y)8(x, y) to the diagonal by 
the very analyticity of the aforementioned expressions. In particular our scheme implies the following 
unique and well-defined definitions of distributions on A4 2 . 


log ra (M 2 q F ) 

a F 


lim /□ 

a—>-0 


lo g n (M 2 (Tp) 


a 


1+a 

F 


n> 0 


(3.32) 


Hereby uniqueness and weak analyticity of /□(log"(Af 2 cr_F)/c^ + “) follow from arguments used through¬ 
out Section 3.3. 


From Proposition 3.5, we know that l/cr)) +Q is weakly meromorphic in a. In order to compute the 
Laurent series, we use the above-mentioned identities for cr and obtain 


1 


n+ 1+a 


2 (n + a) (n — 1 + a) 


(□ + (n + a)f) 


n+a 


in accordance with (3.18) and (3.27). 

Using this and recalling Remark 3.17, we may compute the following Laurent series, where we 
recall that in (3.12) we use the same (arbitrary) constant M present in the logarithmic term of 
(3.28) to correct for the change of dimension and a sufficiently regular function k for later purposes, 


1 1 
(AIk) 2a 

d 11 

1 1 
(. Mh) 2a <j 3 + a 


= 2 (D + /) 


1 log ( M 2 a F ) \ log(fc 2 ) 


aa p 


Uc + /) - 


a F 

2 


(□ + /) □ 

a p Aa p 


a 2 aF 


log 2 (M 2 a F ) 
2(7 F 


□ 


log ( M 2 o f ) + 1 
2(7 F 


io g 2 (fc 2 )(n + f)^ + log (k 2 ) + (□ + /) log ^f F ap ^ 


= ^(□ + 2/)(□ + /)[ — 
8 V aap 


log ( M 2 a F ) \ log(fc 2 ) 


cr F 


(□ + 2 /)(□- 


- - 1 - ((5D + 8/)(□ + /) - 2(D + 2/)/) — + O(a ). 
lo ap 


O(a ), 


(3.33) 


+ 0(a ), 


/) 


1 

a F 
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Note that by means of Lemma 3.18 b) one may explicitly check that the pole terms in these Laurent 
series are local expressions as expected. 

Using the Laurent series, the lowest renormalised powers of ct f may be defined and computed as 4 , 
where we recall once again Remark 3.17. 


log ( M 2 op ) 


lim 

ol— ^0 y 

l 

1 

- pp 

1 1 ^ 

M 2a o 2 F +a 

M 2a a 2 p +a ) 

— lim 

ol — ^0 

(d _ 

1 

1 

d 1 

\da M 2a 

ct 2+ “ 

Op 

PP da M 2, 


Z (J f 


= _!(□ + /) lp g 2 (M 2 v F ) _ D log (M 2 g F ) +1 

4 v ctf 2 <tf 


—) := lim 

oi 


/ 1 


1 


1 


1 


a-s-o V M 2a <ri +c 

r 


- PP 


M 2a a 3 F +a 


- □- 


2(7 p 


(3.34) 


= - ^ (□ + 2/)(□ + /) 1 °g( M2 ^) _ _L ((5D + 8/ )(n + /) - 2(D + 2/)/) — 

8 gf Id gf 


Finally (A|.) ms and (A|>) m _ are defined and computed by expanding the unrenormalised powers A 2 F 
and A'), and replacing the three problematic expressions (3.29) by their renormalised versions (3.34). 


3.5.2. Alternative computation of the renormalised fish and sunset graphs 

As a preparation towards the application of our renormalisation scheme to QFT in cosmological 
spacetimes, we shall now derive an alternative way to compute (A|>) and (A^.) , which is better 

suited for practical computations. We start by stating and proving a few distributional identities. 

Lemma 3.18. The following distributional identities hold. 

a) For any continuous Fq and any twice continuously differentiable F 2 , 

crF 0 S = 0 , o a F 0 5 = 0, F 0 V v *S = ~[F 0 na]5 , 


F 2 DS = [DF 2 ]5 + □ [F 2 }6 - 2X7 a [X7 a F 2 ]S. 


b) 


(□ + /)— =87 T 2 i6 
Op 

(□ + 2 /)(□ + /)—:= lim (□ + 2 /)(□ + f)-^ = 8n 2 i (d - f ) 6 
ctf «-+> off V 3 / 


c) For all n\, n 2 , £ No and n 4 , 715 , n 6 £ {0,1} with n 2 — 77.3 + 774 > — 1, 


log ni (o- F ) (^)%? ( i ) = log n V F ) (* a F r*<T n F 2 


□ log (ctf) = 


□ct — 2 

o F 


V a 


log ns (ct f ) (715 - 7l 6 log” 5 ((7F)) VqCT 


"6 + 1 


d) 


4 Note that we use here a definition of the analytic regularisation of the logarithm in terms of a direct derivative 

rather than a limit of differences like in (3.12). While the two definitions differ up to a constant factor in the principal 
part, they coincide in the constant regular part and thus give the same (<r^ 2 log(M 2 f7ir)) ms . 
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Proof, a) These identities follow from B5 = [B\S for any continuous bitensor B, [a] = 0, [ a a ] = 0 and 
the definition of weak derivatives. 

b) The first identity holds in Minkowski spacetime because 1 /( 87 t 2 ctf) is the Feynman propagator 
of the massless vacuum state. In curved spacetimes (3.30) imply that (□ + f)l/ap vanishes outside of 
the origin and thus must be a sum of derivatives of S distributions. Because o depends smoothly on 
the metric, the coefficients in this sum must be smooth functions of the metric with appropriate mass 
dimension and thus (□ + f)l/ap = cS with a constant c that can be fixed in Minkowski spacetime. 

In order to prove the second identity we recall Remark 3.17 and observe that it is sufficient to 
compute 

t := lim /(□ + /) 

o JK JJ a]+ a 

This expression has for y —> x a scaling degree < 4, vanishes outside of x = y, depends smoothly on the 
metric, is covariant and has mass dimension 6 . Consequently t = cRS where the dimensionless constant 
c can be computed on any spacetime with R ^ 0. Moreover, in any spacetime where / is actually 
continuous in a neighbourhood of the diagonal we have t(x,y) = 8nif(x, x)S(x,y). A spacetime which 
satisfies both properties is (the patch of) de Sitter spacetime defined in conformal coordinates by the 
metric line element 

ds 2 = I (-dr 2 + dx 2 ) 

H z r z 

on (—oo,0) x M 3 , where H is a constant. On this spacetime we have R = 12 H 2 and 


/i 2 := 2 H 2 (t(ti,Xi,t 2 ,X 2 ) = cos 


-l ( T i + t 2 - Oo - 


2 tit 2 


see e.g. [A185], where analytic continuation of cos 1 is understood for time-like separations. From 
this one can infer 

□<r = 1 + 3/iCot(/r) => / = Dg ~ 4 = 6H 2 ^ cot (^) ~ 1 = __ + 0{y 2 ) 

a /j/ 6 

which demonstrates that on de Sitter spacetime / is continuous in a neighbourhood of the diagonal 
with f(x,x) = —R/6. 

c) The distributions on both sides of each equation, considered as distributions in y for fixed x, 
have the same scaling degree < 4 for y —> x and agree outside of the diagonal. Thus they agree also 
on the diagonal as unique extensions. 

d) As in the proof of a) we observe that the potential local correction term on the right hand side 
must be a sum of derivatives of S with coefficients that depend smoothly on the metric because a 
does. Thus the correction term must be of the form cS with a constant c that can be computed in 
Minkowski spacetime. This computation may be performed by using (3.30), the previous statements 
of this lemma, and the following identities which are valid in Minkowski spacetime for any function F 
s.t. F(tjp) is a distribution 


ap\3F(ap) = Dcr fF{cj p) — AF{ap) — 2Vy r7F F(ap ), 


<t pO 2 F(op ) — D“(7 pF^ap ) — 4D F^op') — 4DVvu f F((7f) > 
whereby one finds that c = 0 . □ 

These identities can be used to compute (A 2 F ) and (A'/) in an alternative way under certain 
conditions. 

Proposition 3.19. Let (A A,g) be such that A4 is a normal neighbourhood and let Ap be a distribution 
on M 2 of Feynman-Hadamard form (3.28). Then the following identities hold. 

a) If Ap is a well-defined distribution which is weakly meromorphic in a, then 


(A 2 p) ms = lim 

ot —^0 


M 2 ' 


-Al + “ - 

t 


PP 


M 2 ‘ 


-A 2+c 

l 


flog(87T 2 ) 

167T 2 


5. 
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b) If A p is a well-defined distribution which is weakly meromorphic in a, then 


(Af log (M-Af))„, = lim - '^^6. 

c) If Ap is a well-defined distribution which is weakly meromorphic in a and [v] = 0, then 
(A 3 F ) ms = lim 


«->o l M 2a F 


1 A 3 +“-pp 1 


M 2 ' 


a 3 +cA , i((l + 21 og( 87 r 2 ))i?.+ 1927r 2 [u>]) ^ 
" Af ) + 48(87r 2 ) 2 d ' 


Proof, a) Setting h = 8 it 2 gfAf and k = \J8n 2 /h , we obtain 

1 1 

j 2 + 

7 F 

Using (3.33), [h 2 ] = [ u 2 ] = 1 and Lemma 3.18 a), b) & c) we may compute 


1 A 2+ a _ h- _ 

M 2a F ( 8 tt 2 ) 2 (Mk) 2a <jj +a ' 


lim 


q—> o V M 2a F 


h 2 


1 A 2 +“ - ppiA 2+ “ 


lim 


M 2 ° 
1 1 


1 1 


( 8 tt 2 ) 2 a4o V (Mk) 2a ct 2+ “ PP (Mfc) 2 “ al +a 


h 2 


(8 tt 2 ) 2 \\g 2 , 

\ / \ \ t / ms 


log(fc 2 ) 


Gp 


(D + /) — ) — (A>) ms — 


ilog(87r 2 ) 

167T 2 


S. 


b) In analogy to a), we may compute 


d 1 -A 2+a — nn-_-_ A 2+q 

F PP da M 2a F 


lim i —-—up 

a —>o \ da M 2a 


h 2 ( (log ( M 2 g f ) 


(87T 2 


- tog I w) (k 


i»g 2 («f) 


(□ + /) — 

Gp 


- (APo g (M-A F )), m+ "l i . 

c) This can be proven in analogy to a) and b), whereby one also needs Lemma 3.18 d) and the 
fact that [u] = 0 implies by means of the covariant expansion of bitensors near the diagonal (see e.g. 
[PPV11, Section 5]) that 


v = [w] + ([V a u] - V a [u])cr a + 7 Z v = [V a y]cr a + 7 Z v , 

where the remainder term TZ V vanishes towards the diagonal fast than G a . Thus, the assumption 
[u] = 0 implies that the term in A|, proportional to Gfi 2 log M 2 gf does not need to be renormalised, 
which is crucial for the present proof. The correction term arises from the log/i/(87r 2 ) term in the 
expansion of 

1 1 

(Mk) 2a G 3 F +a 

whose contribution may be computed as 

i (— log^Tr 2 )^ + pit + 87r 2 w;n<7]) S i ((1 + 2 log(87r 2 ))i? + 1927t 2 ['u>]) 

= 8(87t 2 ) 2 = 48(8tt 2 ) 2 ’ 

where again Lemma 3.18 a) & b) prove to be useful. □ 
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3.5.3. A more complicated graph 

In order to show how the proposed renormalisation scheme works for graphs which have more than 
two vertices we discuss the renormalisation of the following triangular graph 

r r := Ap’ |13 Ajr i 23Aj. )1 2 , 

where A p,ij '■= Ap(xi,Xj). In order to apply the forest formula (3.9) to renormalise this graph, we 
note that the forests which correspond to divergent contributions are 

{12}, {123}, {12,123}. 


The renormalisation of rp thus reads 


(Tr)ms — (1 + Rl2 + Rl23 + ^123^12)'^ ^ — (1 + t?123)(l + Rl2)Xp ^ 

In order to illustrate the explicit form of the R, we consider only the most singular contribution to 
Tp a \ namely 


1 


1 


1 


A a ) ._ 

T,0 * _l+o:i3 2(l+ai2) ^-1+0:23 

<7 1 o C'9" 


' 13 O 12 u 23 

where ay, := <7p(x\,Xj). Note that, with obvious notation, (87r 2 )~ 4 w 13 w 2 2 M23ir,o is h 1 fact the only 

7c0 

contribution to -ry which needs to be renormalised. The application of 1 + R \2 to ip 0 has already been 
discussed in the preceding sections and corresponds to the renormalisation of the fish graph. Indeed, 
after setting a 12 , a 23 and 013 to a = aj for I = {1,2, 3} we obtain 


j.(a) 

t T,l 


:= Jim (1 + Ri 2 )tpo = ( ( ) + 0(a ) 


1 


(o-i3) 1+ “ (cr 23 ) 


l+ct ’ 


The distribution (l/<r 2 2 ) ms is a homogeneous distribution of degree S = —A under scaling of X 2 towards 
xi, consequently, tp“j has scaling degree 8 + 4a. 

Owing to Proposition 3.12, we know that tj-j can be decomposed into the sum of a homogeneous 
distribution of degree —8 — 4a and a remainder. Hence, in order to expose the poles of tp“j, we can 
directly apply Proposition 3.9 with m = 1 and Cq = —4a. To this end, we set uq := tpj* and find 


Mi := —4au 0 - E[u 0 = ( ( — 3 - ) + 0(a) 


(vi3) 1+a ( 023 ) 


2+a 


G, 


where G = G(xi, x 2 , 2 : 3 ) is the smooth function introduced in Lemma 3.11. From (3.20) we can infer 
that the principal part of fpj* is 


J-( a ) 

PPH.i = ~ 


4a 


E\ 


G 

V 23 


1 1 
013 023 


whereas the constant regular part can be easily computed as well. Consequently, the renormalised 
distribution 

(*r,o)ms = lim (t { pl - pp 4 Q i) 
can be straightforwardly computed in explicit terms. 


4. Explicit computations in cosmological spacetimes 

The aim of this section is provide pret-a-porter formulae for doing perturbative computations in the 
renormalisation scheme devised in the previous sections for the special case of Friedmann-Lemaitre- 
Robertson-Walker (FLRW) spacetimes. We thus consider spacetimes (A4,g) of the form A4 = I x R 3 C 
M 4 and, in comoving coordinates, 

g = —dt 2 + a(t) 2 dx 2 = a(r ) 2 (—dr 2 + dx 2 ) . 


30 





(4.1) 


Here, t is cosmological time and r is conformal time related to t by dt = adr and 


H := d t log(a) = ^ , 

a z a 


R = 6 (d t H + 2 H 2 ) = ^ . 


We consider here the spatially flat FLRW case for simplicity. Note that these spacetimes are normal 
neighbourhoods so that (3.28) can be considered as a global expression and all Feynman amplitudes 
can be analytically regularised without the need of introducing partitions of unity such as in Remark 
3.1. 


4.1. Propagators in Fourier space 

In comoving coordinates with conformal time, the Klein-Gordon operator reads 

P = -□ + + TO 2 = ^<9 2 - V 2 + ^ Ra 2 + m 2 a 2 ^ a(r). 

It is convenient to employ Fourier transformations with respect to the spatial coordinates in order to 
expand quantities in QFT on FLRW spacetimes in terms of mode solutions of the free Klein-Gordon 
equation 




XfcO) 


Akx 


(27t) 2 a(r) 


where the temporal modes Xk( T ) satisfy 


d 2 + k 2 + m 2 a 2 + ( £ - ^ ) Ra 2 ) Xk(r) = 0 


(4.2) 

(4.3) 


and the normalisation condition 

XkdrXk - XkdrXk = * • 

Here, k \= |fc| and 7 denotes complex conjugation. 

In particular, we can use the mode expansion in order to give explicit expressions for the various 
propagators of the free Klein-Gordon quantum field in a pure, Gaussian, homogeneous and isotropic 
state (see [LR90, PilO, Zsl3] for associated technical conditions on the mode functions). To this 
avail, we define 


1 


A|j(xi, x 2 ) =: lim / d A k A tt (Ti, t 2 , k) e' 

40 87r 3 a(n)a(r 2 ) J R3 


ik{x i —X 2 )—ek 


(4.4) 


where Aj stands for either A + (two-point function), Ar/a (retarded/advanced propagator) or Aj? 
(Feynman propagator). See Section A.l for our conventions for these propagators and their relations. 
Recall that our renormalisation scheme preserves invariance under spacetime isometries and thus we 
know that renormalised powers of the Feynman propagator may also be written in the form (4.4). 
The Fourier versions of the single propagators read 


A+(ti,t 2 , k) = Xfe(ri)x fc (T 2 ) , A_(n,T 2 ,fc) = A + (n,r 2 , k ), 

A F {ri,r 2 ,k) = 0(ti - t 2 )A + (ti,t 2 , k) + 0(r 2 - ri)A_(n,r 2 , k ), 
A R /A(n,T 2 ,k) = =F*©(±(ri - r 2 )) (a + (t 1; t 2 , k) - A_(n, r 2 , k)') 


(4.5) 


whereas by the convolution theorem, we have the following Fourier versions of products and convolutions 
of multiple propagators, provided those products and convolutions are well-defined. Defining 


[A#! *4 Aj 2 ] (x, y) ■.= / d 4 x^g A h (x 1 ,x)A h (x,x 2 ) 
J M 


A#! *1 A(J 
A#! *3 Aj 


(n ,r 2 ,k) := J dr a(r) 2 A h (n,T,k) A h (r,T 2 ,k) 

(ri,r 2 ,fc):= / d 3 p A h (n, r 2 ,p)A|j 2 (■n, r 2 ,\k - p\) 
J R3 V ' 


(4.6) 
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we have 


IW Ti ^ fc ) = 


1 


\n— 1 


((27r) 3 a(ri) 2 a(r 2 ) 2 ) 

Ajjj *4 • • • *4 = Ajj *i • • • *i Aj n . 


A #1 *3 • • • *3 A|(„ (n,r 2 ,fc) 


(4.7) 


Choosing a pure, Gaussian, homogeneous and isotropic state G of the quantized free Ivlein-Gordon 
field on a spatially flat FLRW spacetimes amounts to choosing a solution of (4.2) and (4.3) for each k. 
In order for G to be a Hadamard state the temporal modes Xk have to satisfy certain conditions in the 
limit of large k which are difficult to formulate precisely. Heuristically, a necessary but not sufficient 
condition is that the dominant part of Xk for large fc, when the mass and curvature terms in (4.2) are 
dominated by k 2 , is ^=e _zfeT , i.e. a positive frequency solution. Note that the retarded and advanced 

propagators are state-independent and thus A r/a(tii t 2 i k) is independent of the particular Xk chosen 
for each k. 

4.2. The renormalised fish and sunset graphs in Fourier space 

In perturbative calculations at low orders we encounter (pointwise) powers of A± and A F . While the 
powers of A± are well-defined if G is a Hadamard state on account of the wave front set properties 
of these distributions, we need to renormalise the powers of A F by means of the scheme developed 
in the previous sections. In order to be useful for explicit computations in FLRW spacetimes, we 
have to develop a spatial Fourier-space version of this scheme. Having in mind the application to 
</> 4 theory, we shall compute (Ap’)g ls (n, t 2 , k) for n = 2, 3. The difficulty in achieving this is that, to 
our knowledge, despite of the large symmetry of flat FLRW spacetimes, neither a nor the Hadamard 
coefficients u, v and w may written in a tractable form which can be Fourier-transformed easily. Our 
strategy to circumvent this problem is the following. 

Computational strategy 

a) For a general mass m and coupling to the scalar curvature £ and a general homogeneous and 
isotropic, pure and Gaussian Hadamard state G, split A F as 

Af = A Ft0 + d , d := A F - Af )0 , (4.8) 

where A F q must satisfy the following conditions. 

• Ap t o is explicitly known in position space and Fourier space. 

• A F q is of the form 

Ago = g ~2 (^~ + lo S ( M 2 a F )^j + w 0 , 

with ito = u, i.e. it agrees with A F in the most singular term but not necessarily in the 
subleading singularities. 

• [t>o] = 0 and A p 0 is weakly meromorphic in a such that (A|, i0 ) , (A|. 0 log (M -2 A^o)) 

and (A|. 0 ) ms may be computed with Proposition 3.19. This is crucial for preserving the 
explicit knowledge of Ago in position space in the renormalisation procedure, so that one 
may hope to compute the Fourier transforms of the renormalised powers. 

b) With these assumptions on Ago it follows that the renormalised fish and sunset graphs may be 
computed as 

(Af)L = (Ago)ms + 2A F' 0 d + d 2 (4.9) 

(Af)L = (Aigo)ms + ^ (A>,o^) ms + 3A F fid 2 + d 3 

because the non-renormalised terms in the above formulae are distributions with scaling degree 
< 4 for y —> x and thus can be directly and uniquely extended to the diagonal. 
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c) (A|i )0 ) and (A f>0 ) may be computed with Proposition 3.19 as anticipated. In order to 
compute (A|, 0 d) , we further split d as 


d — d\ o?2 , d\ 


M log (m 2 A Fi o) 
87T 2 


d 2 := d — d\. 


(4.10) 


Because v = [u] + 0(a a ), d\ contains the leading logarithmic singularity in d (and thus A F ) 
which is the only logarithmic singularity relevant for the renormalisation of the sunset graph. 
Consequently 


(A F,o d ) ms - ( A f,o 1o S i M 2 AF.o)) ms + d 2 (A po) 


and thus Proposition 3.19 can be applied again. 


(4.11) 


d) Due to the symmetry of FLRW spacetimes and the assumption that the pure and Gaussian 
Hadamard state Q is invariant under this symmetry, [u] and [ic] do not depend on the spatial co¬ 
ordinates. Given that one succeeds to compute the spatial Fourier transforms of log ( M~ 2 Af,o ), 

( A F,o) ms > ( A F,o lo g ( M_ 2 A Co)) ms and (A 3 F0 ) (A F )™ s (n, r 2 , fc) for n £ {2,3} may be com¬ 

puted by means of the convolution identities (4.7), since the Fourier transforms of A F0 , d\ and 
d 2 are known by construction. 


In order to follow the computational strategy outlined above, we first compute [a] and [in]. Indeed, 
the coinciding point limit of the Hadamard coefficient v reads (see e.g. [HalO, Section III.1.2] for 
details) 


v = 


77,2 + (£ - |) R 


(4.12) 


Moreover, using the method of [SclO] to compute a spatial Fourier representation of the Hadamard 
parametrix Hp - here considered as (3.28) with w = 0 - in FLRW spacetimes, one can compute (see 
the review in [Del3] and a related method in [PilO] for the conformally coupled case) 


[in] = lim (Ap(x,y) — Hp(x,y)) 


x^y 


(2 F ) 3 a 2 J d3k lxfc(r)l " 


k 2 + a 2 m 2 + a 2 (£ — |) R 


(4.13) 


1 ^ + (i-i)fl)(2 7 -i + i„ g ( m2 + ( «G)«' 


16-7T 2 


R 


2 M 2 


36(8tt 2 ) 


21 ’ 


where 7 is the Euler-Mascheroni constant and Hp is taken with the mass scale M inside of the 
logarithm of a 5 * * . 

As anticipated we see that [n] and [in] are functions of time only (recall (4.1)). Moreover, we see 
that [n] = 0 for a conformally coupled (£ = f) massless scalar field. Thus, in order to pursue our 
computational strategy, we should look for a candidate for A F o among the Feynman propagators 
in suitable states of this theory. In fact, choosing the conformal vacuum state of the massless 
conformally coupled scalar field does the job. The conformal vacuum is given by choosing the modes 
Xkij) — e~ lkT /V2k, and thus the Feynman propagator A F q in this state is of the form 


a f, 0(^15 ^ 2 ) = 


87r 2 a(ri)a(r 2 ) 0 \f,m(^i, z 2 ) 


A F,o( r li r 2, k) = 


a -ik\r 1 -T 2 \ 

2k 


(4.14) 


Here, and in the following, the index m indicates quantities in Minkowski spacetime, in particular 
ctm(£i,;E 2 ) = 5(^1 — X 2) 2 ~ 5(17 — r 2 ) 2 . A fo is weakly meromorphic in a because the massless 

5 Note that one may take instead of the function F(k) = 1/(2^/ k 2 + a 2 m 2 + a 2 (£ — jl) R) in (4.13) any distribution 

F'(k) such that F'(k) — F(k) is 0(k~ 5 ) for large k and integrable. By taking e.g. F’(k) = 1/(2 k) — 0(fc — am)(a 2 m 2 + 

a 2 — {) _R)/(4fc 3 ) one may cancel the log R term outside of the integral. 
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vacuum Feynman propagator in Minkowski spacetime has this property and the conformal rescaling 
by a does not violate it. Thus, we may follow our computational strategy and compute (Ap 0 ) , 

(Ap 0 log (M~ 2 Af,o)) and (A|. )0 ) by means of Proposition 3.19. This is easily done using (3.33) 
for <jf,m rather than of and h = ^/ 87 r 2 a(ri)a(r 2 ) = \/8tt 2 ci <S> a. The results are 


( a f,o)L = iim 


a—>-0 \ 


— (A F , 0 ) 2 +“-pp^(A F ,o) 2+ “ 


ilog(87r 2 ) 

167T 2 


8 


_ lim 1 ( 1 1 1 1 \ 
a->o ( 87 r 2 ) 2 a 2 ® a 2 ( (Mv / 87 r 2 a <g> a) 2 “ (M\/8n 2 a ® a) 2a o J 

1 + 21og(a) 1 \og(M 2 o F ,m) 

1 (2 2 4 o/o 2\2 2 xrx 2—M ’ 

107r^a 4 2(87 v z ) z a z (g) a z &f,m 


i log(87r 2 ) 
167T 2 


s 


(A| 0 log(M 2 A F , 0 )) ms = lim 




ms a-s-o \ da M 2a 

2 + 21 og(a 2 87 r 2 ) + log 2 (a 2 ) . 


i8t, 


1 


+ 


32-7r 2 a 4 4 ( 87 r 2 ) 2 a 2 ® a 2 

1 + log( 87 r 2 )a ® a log (M 2 o f ,u) 

-Ur- 


□f 


flog (87r 2 ) 

■ 32 jt 2 S 

log 2 (M 2 <tf,r 
i- 

ff F,M 


2 (87r 2 ) 2 a 2 (g> a 2 




and 

( A F.o)ms = lim n 


“Ao V M 2a ( Af ’°) 3+ pp M 2 “ ( Af ’°^ 3+ ) + 


i ((1 + 21og(87r 2 ))f? + 1927 t 2 [w]) 


(15 + 12 log (a)) D m + 6 (D M log(a)) + 2 (d 2 a)/a . 


48 (87r 2 ) 2 a 6 


7(5m — 


48(8tt 2 ) 2 


8(87r 2 ) 3 a 3 <g> a 3 


□i 


log (M 2 o f ^ 

i “ 

&F, M 


where we have used 8 = Sm/cl 4 , fu = 0 and the fact that, by (4.13), 87 t 2 [ujo] = —R/ 36 for the conformal 
vacuum state of the massless, conformally coupled scalar field. 

Using these results as well as the Fourier representation of 1/of,u (4-14) and log (M 2 of, m) (A.2), 
and convolution identities, we can finally obtain the Fourier versions of the renormalised powers of 
A p,o- For instance, we find for (A|, 0 ) 


( A f,o) (n,T 2 ,k) 

\ / ms 


l + 21 og(a(ri)) 

- 6(n - r 2 )- 


1 


- (d 2 

167r 3 a(ri)a(r2) Tl 


k 2 ) [ d 3 P tU±r 


16a(ri) 2 7r^ 

, *l r l ^21 

ren 


(4.15) 


1 


2 p 2 


2 | k-p\ 


a -i(p+\k-p\)\Ti-T 2 \ 


where the appearing renormalisation of 1/p 3 is defined in (A.l). Note that the p-integral has no conver¬ 
gence problems for large p because one may write the potentially dangerous — i\r\ — T 2 \e~ 2lp ' Tl ~ T2 '/p 
contribution as d p (e~ 2lp ' Tl ~ T2 ' / (2p 2 )) plus an 0(p~ 3 ) term. Regarding the convergence for small p we 
observe that the integral is manifestly convergent if A; 0, thus yielding a well-defined distribution in 
k on R 3 \ {0}. The scaling degree of this distribution is easily seen to be 1 < 3 and thus a unique 
extension towards the origin exists. In practical terms this means that the integral for k = 0 may be 
computed as a limit k — > 0 of the integral with nonvanishing k without any renormalisation. 


4.3. Example: the two-point function for a quartic potential up to second order 

In order to compute the analytic expressions corresponding to the graphs in Figure 2, we may use 
the Fourier versions of the appearing propagators (4.5), (4.15), and the analogous expressions for 

( A|. 0 log (M~ 2 A|, 0 ) 1 (t - !, 7"2 , /c) and (A|, 0 ) (r-j,r 2 , k) the explicit form of p(x) = 3A w(x,x) in 
(4.13), as well as (4.9), (4.11) and the identities for products and convolutions (4.6), (4.7). Note that 
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^{x) is in fact only time-dependent because was chosen homogeneous and isotropic. Thus the 
integrals with /z-vertices can be computed partly with the above-mentioned identities by means of 

(1 <g>/i)Ajt(Ti,T 2 ,/c) = /Lt(r 2 )Att(ri,T 2 ,fc), (/x <S» 1)A|j(ti, r 2 , k) = Aj|(ri, r 2 , k). 


Similarly, the bubbles in the third line of Figure 2 contribute only time-dependent vertex factors which 
can be computed as 

fyj( T ) := f d,Tid 3 Xi a(r 1 ) 4 /r(r 1 )Ajj(r, ri, * — £ 1 ) = — [ dr x a(ri) 3 ^(n)A|j(r, r 1; 0) 

Jm a \ T ) Ji 

where A# is either A+ or (A|>) ms . 

With these preparations, we can compute e.g. the first graphs of the fourth and fifth line in Figure 
2 in Fourier space as 


A^ *4 {{hp <8> 1)A + ) = Ar *i ({hp <g> 1)A+)) 


= / dT 3 dT 4 a(r3)a(T4) 3 /r(T4)A7 ? (Ti,T3,fc)A + (T3,r 2 ,fc)(A|,) ms (r3,r4,0) 
Ji 2 


and 


A_R *4 (Aj?)^ *4 A + = A R *1 (A F )ms *1 A + 


= / dr 3 dr 4 a(T 3 ) 2 a(r 4 ) 2 A i? (Ti, r 3 , fc)(A|) ms (r 3 , r 4 , fc)A + (r 4 , r 2 , fc). 
Ji 2 


4.4. More complicated graphs on cosmological spacetimes 

In order to compute the Fourier transforms of more complicated graphs on FLRW spacetimes, one 
can use a strategy generalising the one employed in Section 4.2. Namely, one again decomposes 
the Feynman propagator A p into several pieces which capture the relevant singularities and can be 
expressed in terms of the conformal vacuum Feynman propagator A^ 0 whose explicit form in position 
and Fourier space is well-known in contrast to the form of a itself. The corresponding decomposition 
of general Feynman amplitudes it is straightforward. The only non-trivial step is to generalise 
Proposition 3.19 to the case of general amplitudes, i.e. to compute the difference between the minimal 
subtraction scheme used in conjunction with either analytically regularising powers of a directly or 
analytically regularising powers of the full propagator A p 0 . However, we do not foresee any problems 
in obtaining such a generalisation by proving versions of Lemma 3.11 and Proposition 3.12 for A^o 
rather than a. 

In fact, one can also skip this last step by taking a rather pragmatic approach and working directly 
with the renormalisation scheme consisting of decomposition in Aj?o, analytic regularisation of powers 
of this propagator and minimal subtraction of the principal parts. This scheme, clearly applicable only 
to conformally flat spacetimes, satisfies all properties proved in Proposition 3.14, with two exceptions. 
It is not obvious whether the Principle of Perturbative agreement with respect to generalised mass 
perturbations holds for this scheme, whereas locality and covariance of course only hold in the sense 
restricted to conformally flat spacetimes. In this respect it is essential that the Feynman propagator 
of the conformal vacuum Ap o on conformally flat spacetimes is manifestly “geometric”, because the 
corresponding propagator of the massless Minkowski vacuum has this property. 

5. Summary and outlook 

In this work, we have introduced a renormalisation scheme on curved spacetimes consisting of a 
particular analytic regularisation of the Feynman propagator, and thus of all Feynman diagrams, and 
a minimal subtraction of the principal (pole) part of the resulting meromorphic expressions. We have 
argued that this scheme has all properties that a physically meaningful renormalisation scheme on 
curved spacetimes should have and that it is in fact a particular form of differential renormalisation. 
The renormalisation scheme discussed in this work has the advantage that it is 
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a) directly applicable to spacetimes with Lorentzian signature, 

b) manifestly (local and) covariant, 

c) manifestly invariant under any spacetime isometries present, 

cl) capturing correctly the non-geometric and non-unique state-dependent contribution of Feynman 
amplitudes and not only the geometric divergent part, which is unique up to finite renormalisa¬ 
tions, 

e) well-suited for practical computations, e.g. in cosmological spacetimes, 

f) constructed to all orders in perturbation theory, 

g) and mathematically rigorous. 

To the best of our knowledge, other renormalisation schemes on curved spacetimes discussed in 
the literature such as dimensional regularisation, local momentum space methods, zeta-function 
regularisation, heat-kernel techniques, generic Epstein-Glaser renormalisation and, on cosmological 
spacetimes, dimensional regularisation only with respect to spatial variables, lack at least one of the 
features listed above. 

In order to demonstrate the practical applicability of the scheme, we have computed several 
examples on generic curved spacetimes. Moreover, we have shown how explicit computations in 
cosmological spacetimes can be done, in particular, how the renormalisation scheme initially defined in 
position space can be interpreted in terms of quantities Fourier transformed with respect to comoving 
spatial coordinates. 

We have discussed the renormalisation scheme only for scalar fields in four spacetime dimensions, 
however, the extension to other spacetime dimensions is straightforward. Moreover, as the analytic 
regularisation discussed in this work consists of regularising only inverse powers of the squared geodesic 
distance, it can be straightforwardly generalised to field theories with higher spin, with and without 
gauge-invariance. In particular, spinorial quantities can be directly regularised without the need to 
worry about their dependence on the dimension such as in dimensional regularisation. Finally, we 
expect that a generalisation of the scheme introduced in this work to gauge theories yields a scheme 
which preserves or can be modified to preserve the local gauge symmetry in theories which are free of 
anomalies. 


Acknowledgments. The authors would like to thank Klaus Fredenhagen and Markus Frob for 
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Forschungsgemeinschaft (DFG). 

A. Conventions and computational details 

A.l. Propagators of the free Klein Gordon field and their relations 

A+(x,y) = (<j>(x)<j>(y))n = (<j>(x) * <j>(y)) n = {<i>(y)<j>(x ))n = A -(x,y) 

A(x,y) = - (A +(x,y) - A -(x,y)) = A R (x,y) - A A {x,y) 

l 

A R/A = ±©(±(G - tyj) A(x, y) 

A F (x,y) = (T (<l>(x)<t>(y)))n = {<f>(x) - T 4>{y))n = ©(G - t y )A+(x, y) + Q(t v - t x )A_(x,y) 

=> A F (x,y) = ^ (A+(x,y) + A_(x,y)) + l - (A R (x,y) + A A (x,y)) 

= A + (x,y) + iA A (x,y) 

= A_(x,y) +iA R (x,y) 
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A.2. Fourier transform of the logarithmic term on Minkowski spacetime 

In order to compute the Fourier-transform of log (M 2 <tf,m), we recall that the Feynman-propagator of 
the Klein-Gordon field with mass m in the Minkowski vacuum is given by 


— hill 


to (27 t) 3 


d 3 k 


= —' iVk 2 +m 2 \T 1 -T 2 \ „ik(x—y) p -ek 


2 VP 




1 / 2 m 2 

= ^—2 A / —- K 1 ( \/2m 2 a FM 

a F,M 


8ir 2 

1 

8n 2 


1 m 2 / m 2 a F M \ , / e 2 ^m 2 a FM \ m 2 f 5m 2 a FM \\ 4 

+ ^ v + 4 J l0 ® i 2 + ~8 jj +0(m) ' 


where K\ is a modified Bessel function and 7 is the Euler-Mascheroni constant. Using this, we find 
log (M 2 ct f ,m) 


= lim ( I 67 r 
m—¥ 0 
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dm 2 


-log 


2 M 2 


= — lim lim — 

m -¥0 l etO 7T , 

= — lim — 

<40 7r 


d 3 t l + iy / k2 + m2 \ — ~ T2 1 '■“ i 
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<40 7T 
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d 3 k 
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ym-^O 

, *|n - 72 1 

\ 

2 k 2 

d 3 k 

6G 


2 {k 2 + m 2 ) 2 
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e -i v / fc 2 +m 2 |T 1 -r 2 | e ifc(x-y) g-efc _|_ ^ 0 g ^ 
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2 U 3 


flog(n - r 2 , fc) := -V8tt 
is the sought-for spatial Fourier transform of log (M 2 a F} m)- 
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